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The growing interest in the study of the chemical composition of cellular 
formed elements has led to wide-spread research on techniques for demon- 
strating and estimating, in situ, the presence and amounts of the various 
components found by macrochemical analysis. In short, this would in- 
volve the use of a.quantitive chemical reaction and a means of accurately 
measuring the end-product of the reaction. Such techniques for measuring 


nucleoproteins of cell nuclei were recently described by Pollister and Ris.! 
Using these same methods, a study of the nature of the Feulgen nucleal re- 
action was undertaken. 


The Feulgen reaction may be considered as occurring in two steps: first, 
the splitting of the linkage between the purine bases and the desoxypentose 
of the nucleic acid by means of mild acid hydrolysis, and second, a chemical 
reaction between the aldehyde groups, thus formed on the sugar, and the 
Schiff reagent (decolorized fuchsin) resulting in the synthesis of a new dye 
compound. It was noticed by Bauer,” and confirmed by Hillary,* that the 
intensity of color produced in chromatin of cell nuclei is related to hydroly- 
sis in such a way that as time of hydrolysis is increased there is a more and 
more intense staining until a maximum point is reached; while beyond this 
maximum, longer hydrolysis results in a decreased staining. The in- 
creased staining could easily be explained by the progressive splitting off of 
bases, as a consequence of which more and more sugar aldehydes are made 
available for the reaction, which would lead one to expect a maximum when 
all the purine-desoxyribose linkages have-been broken. But the decrease 
in color after the maximum intensity of stain is reached is quite unexpected. 
Bauer offered two possible explanations for the latter phenomenon: first, it 
might be due to a continued decomposition of that part of the nucleic acid 
molecule remaining after the purine bases have been split off; or second, it 
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might be due to the destruction of the protein-nucleic acid linkage allowing 
the latter to diffuse from the chromatin. Hillary subsequently showed 
that the decrease in intensity could not be due to a chemical degradation of 
the aldehyde or color-producing groups in the sugar, because a solution of 
thymus nucleic acid in a test tube is still capable of reaction to produce color 
after eight hours of hydrolysis. Furthermore, Hillary hydrolyzed nucleic 
acid agar blocks from five to thirty minutes, and after testing both block 
and hydrolyzing medium with the Schiff reagent found that the staining 
capacity of the latter was roughly the reciprocal of that in the agar blocks. 
This last experiment suggests that one result of prolonged hydrolysis is to 
change the nucleic acid into a more readily diffusible form. - 

In order to reach a clearer understanding of what might be 
happening to the nucleic acid molecule as time of hydrolysis is increased, 
Professor Arthur W. Pollister suggested that an attempt be made to follow 
the fate of the molecule by separate photometric determination of each of 
the three constituents of the individual mononucleotides. Cartilage cells 
from the head region of frog tadpoles were selected as-especially suitable for 
these measurements, since in these cells the nuclei are nearly perfect spheres 
and the cytoplasm contains little ultra-violet absorbing material. The di- 
ameter of these nuclei is about seven microns. The cartilage was fixed in 
Carnoy’s acetic alcohol (1:3), and sections were cut at ten microns. Only 
whole nuclei were measured, and corrections were made for the cytoplasmic 
absorption and for the non-specific light loss as determined in appropriate 
blanks. The fate of the nucleic acid bases was followed by means of ultra- 
violet absorption, that of the desoxypentose by measuring the amount of 
dye produced by the Feulgen reaction performed in the manner described 
by Stowell,‘ and that of the phosphoric acid by overstaining with methyl 
green and removing most of the dye not chemically combined, by prolonged 
washing in absqlute isopropyl alcohol. 

Bauer constructed curves based only upon visual estimates of Feulgen 
dye intensity; and the similar curve of Fig. 1, which is based upon actual 
photometric measurements of absorption by the dye in nuclei, agrees fully 
with his semi-quantitative observations. There is a peak at 12 minutes 
which is followed by a decline until at 24 minutes the extinction value is 
barely measurable. 

The change in concentration of purine and pyrimidine bases throughout 
hydrolysis is shown by the curve based on the ultra-violet absorption. 
Feulgen® found that mild acid hydrolysis, like that in the nucleal reaction, 
readily breaks the purine-sugar linkage, while by contrast, the pyrimidine- 
sugar linkage is broken only after much more drastic treatment. It has 
been recently shown that most, perhaps all, desoxypentose nucleic acids 
contain equal amounts of purines and pyrimidines (Mirsky and Pollister‘) ; 
consequently, in the Feulgen reaction the time when the ultra-violet ex- 
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tinction has reached a value one-half that of the unhydrolyzed nuclei would 
be expected to correspond to that point at which all the purines have been re- 
moved by hydrolysis. The extinction for unhydrolyzed nuclei was found to 
be 0.140. A measured value of 0.070, one-half that of the unhydrolyzed 
nuclei, was found at 12 minutes hydrolysis. As we saw above, this was the 
moment when the amount of regenerated fuchsin dye was at its maximum. 
It appears likely therefore, that for the first 12 minutes’ hydrolysis the 
Feulgen reaction is proceeding in the accepted manner, i.e., by the progres- 
sive splitting off of purines producing in the desoxypentose more and more 
free aldehyde groups which become available for restoring color to the de- 
colorized fuchsin. 

Turning to a consideration of the effect of longer hydrolysis, it can be 
seen from the ultra-violet absorption curve that beyond 12 minutes the 
amount of purine and pyrimidine decreases steadily. Finally after 24 
minutes’ hydrolysis, at the time when we have seen that the Feulgen re- 
generated dye indicated practically no reacting sugar, the bases likewise 
are no longer present. ; 

The original nucleic acid of the nuclei contained, presumably, one phos- 
phoric acid radical for each desoxypentose and base residue. The methyl 
green curve follows the history of this component during the course of the 
standard hydrochloric acid treatment. At the end of 12 minutes’ hydroly- 
sis the nuclei contain the full amount of phosphoric acid. Thus, during a 
period when up to one-half the base content is being removed from the 
nucleic acid, and while the Feulgen reaction with the desoxypentose is con- 
sequently increasing to a maximum, the third component of the nucleic 
acid has remained unchanged. This is strong evidence that, during a hydrol- 
ysis which is optimal for production of the Feulgen nucleal color, the change in 
the nucleic acid molecule involves only the splitting off of purine bases, as orig- 
inally postulated by Feulgen. The phosphoric acid curve likewise sheds 
most interesting light on the effect of longer hydrolysis. Beginning at 12 
minutes the phosphoric acid decreases rapidly and at 24 minutes the amount 
of methyl green bound by the nuclei is barely measurable. Taken with the 
evidence from fuchsin and ultra-violet absorption, this shows conclusively 
that the events from 12 minutes on are very different from those during the 
earlier part of the hydrolysis; for the overhydrolysis leads to the destruction or 
loss from the nuclei of all three components of the nucleic acid molecule. Inthe 
light of Hillary’s experiments (see above) it is probable that the latter is the 
correct explanation; the overhydrolysis has altered the nucleic acid so that 
it freely diffuses out of the chromatin. It would seem that this very likely 
involves some depolymerization. Another process that may be expected to 
be taking place is a breaking of the nucleic acid-histone salt linkage; for 
Carnoy’s reagent has been found by Pollister and Ris to be unique among 
cytological fixin~ fluids in that it preserves the nuclear histones in such a 
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form that they are readily removable by chemical means (e.g., by solution 
in the divalent mercury-sulfuric acid Millon reagent). 

This brings us to the question of the relation of fixation to the Feulgen 
reaction, which was exhaustively studied by Bauer. He found sharp- 
peaked nucleal reaction curves, like that above, after fixation with such 
fluids as Carnoy and acetic sublimate. A similar curve has also been found 
for a trichloroacetic sublimate fixing fluid (Pollister unpublished). By 
contrast, Bauer found that when tissue is fixed in fluids containing chromic 
acid (such as Allen’s Bouin, San Felice, Flemming, Flemming-Heitz, and 
Champy) instead of a sharply restricted moment of optimal hydrolysis 
there is a long optimal period during which there is no change in stain in- 
tensity; although very prolonged hydrolysis eventually does lead to a drop 
in stain intensity. It is interesting that the fixing fluids which Bauer found 
to give sharp-peaked curves upon hydrolysis are those which fix tissue in 
such a way that all, or nearly all, the histone is readily removable; while 
from chromic acid fixed material the same treatment removes little or none 
of the protein (Pollister unpublished). Experiments are iow being per- 
formed to measure the amount of histone at different stages of hydrolysis. 

A compelling reason for undertaking the present study was to examine 
the possibility of using the Feulgen reaction for quantitative determination 
of desoxypentose nucleic acid in a cytological preparation. The following 
considerations indicate that, under proper conditions, the intensity of the 
Feulgen reaction is an accurate measure of the amount of desoxypentose 
nucleic acid in a nucleus. The stain intensity at the peak is fairly ac- 
curately reproducible for one type of tissue, and the extinction value at the 
peak after chromic acid fixation seems to be close to that after Carnoy’s 
fluid. As we have seen above, this peak probably represents a moment 
when all the purines have been removed. From the extinction value of the 
optimal Feulgen reaction, one can readily compute the amount of dye in an 
average nucleus (standard value for regenerated fuchsin, 4.4 X 10~-® 
gram per cc., 1 cm. path, extinction at 530 my, 0.745). It is Wieland and 
Scheuing’s’ view of the nature of the Schiff reaction, that two parts of each 
decolorized fuchsin molecule react with two aldehyde groups to restore the 
color. Hence, if it is assumed that the color in the nucleus after optimum 
hydrolysis has resulted from the reaction of one dye molecule with two de- 
soxypentose residues to which a purine was originally attached, the total 
amount of nucleic acid in an average nucleus can be computed. A wholly in- 
dependent computation of nucleic acid in a cartilage cell nucleus may be 
made from the extinction, near the ultra-violet absorption peak, of nucleic 
acid in an unhydrolyzed preparation. By calculation from ultra-violet 
measurements, the average cartilage nucleus contains 2.44 X 10~-!? gram of 
nucleic acid; from the Feulgen reaction the average nucleus contains 2.39 X 
10-”* gram of desoxypentose nucleic acid. A certain amount of justi- 
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fication for the use of methyl green staining as a measure of nucleic acid 
phosphoric acid comes from the fact that if one assumes that each of the di- | 
basic methyl green molecules unites with the phosphoric acid of two nucleo- 
tides the amount of nucleic acid in the average cartilage nucleus is computed 
to be 2.44 X 10-'* gram. 

It should be added that the above figures are of the same order of magni- 
tude as previous determinations of total desoxypentose nucleic acid in a 
single nucleus. Pollister and Ris report that the nucleus of the calf thymus 
lymphocyte contains 1.0 X 10—- gram, from cytological determination by 
ultra-violet absorption; while calculation from macrochemical data on the 
lymphocyte nucleus gave practically the same figure, 1.1 X 10~—!* gram. 
Caspersson® found 11.2 X 10~—!* gram of nucleic acid in the large grass- 
hopper leptotene spermatocyte nucleus, by the ultra-violet cytological 
method. From the macrochemical data of Zittel and O'Dell’ the amount 
- of desoxypentose nucleic acid in the nucleus of a single bull spermatozoon 
can be computed to be 4.6 X 10-1? gram. 

Summary.—By photometric measurements the base, sugar and phos- 
phoric acid of the desoxypentose nucleic acid of cartilage nuclei have been 
followed throughout the standard hydrolysis of the Feulgen nucleal re- 
action. The results indicate that at optimum hydrolysis time one-half the 
total base content has been removed; presumably this represents the pur- 
ines. During post-optimal hydrolysis the whole nucleic acid molecule dis- 
appears from the nuclei. 

The possibility of using the Feulgen reaction for quantitative determi- 
nation of desoxypentose nucleic acid is briefly considered. 
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DIETARY FACTORS IN THE UTILIZATION OF HOMOCYSTINE* 


By G. W. Kipper AND VIRGINIA C. DEWEY 
BIOLOGICAL LABORATORY) AMHERST COLLEGE, AMHERST, MASSACHUSETTS 


Communicated by E. W. Sinnott, January 15, 1948 


Introduction.—It has been known for some time, mainly ‘due to the in- 
vestigations of du Vigneaud, ef al., !.? that homocystine can replace 
methionine in the diet of the rat provided a methyl donor, such as choline 
or betaine, is present. The recent work of McKittrick* has pointed to the 
same conclusion regarding the chick. 

Evidence for the ability of the rat to synthesize methionine from homo- 
cystine in the absence of choline, betaine or any other known methyl donor 
was first presented by Toennies, et a/.4 They found that rats, maintained 
on a rich pre-experimental diet and then on a depletion diet, lacking cho- 
line and sulfur containing amino acids, responded to the addition of choline 
plus homocystine or to homocystine alone. They discuss the influence of 
intestinal microérganisms on their results and attempt to minimize this in- 
fluence by the administration of sulfasuxadine*®. Even after prolonged sul- 
fonamide treatment their rats showed a positive growth response to homo- 
cystine but only in the presence of a liver fraction. This liver fraction was 
essentially free of methionine and choline. The authors concluded that 
liver contained a factor or factors making possible the methylation of 
homocystine from sources other than choline or betaine. 

In the course of our investigations of the biochemistry of the ciliated pro- 
tozoan, Tetrahymena, we found that homocystine could replace methionine 
in its diet. Choline was present as was a crude growth factor preparation 
(Factor II) from liver. When Factor II was prepared from plant sources, 
homocystine alone did not support growth. In addition, subsequent ex- 
periments have shown that exogenous choline does not influence the 
methylation of homocystine but that our Factor II preparation from liver 
contained the agent described by Toennies and Bennett, while Factor II 
prepared from plant material (Cerophyl) was devoid of the Toennies and 
Bennett factor. 

The purpose of this report is to present our findings, which support those 
of Toennies and Bennett, on the occurrence of a factor or factors in liver 
which allow the methylation of homocystiae to methionine in the absence 
of exogenous choline and betaine. The possibility of unknown microér- 
ganisms vitiating the results, always present in work done with verte- 
brates, is obviated in our work. Tetrahymena is grown in pure (bacteria- 
free) culture and our results reflect the animal’s metabolism only. More- 
over, it should be remembered that Tetrahymena possesses a rather typical 





82 BIOCHEMISTRY: KIDDER AND DEWEY Proc, N. A. S. 


animal biochemical pattern so that comparisons of our results with those 
on other animals are permissible. 

Experimental.—The organism used in this investiagtion was the ciliated 
protozoan, Tetrahymena geleii W grown in pure (bacteria-free) culture. 

In our previous studies on the biochemistry of this organism 7—'? growth 
responses were measured by counting the number of ciliates, after appropri- 
ate dilutions, in measured volumes of media.. The results were expressed 
as organisms.per ml. of fluid. Various attempts to measure growth tur- 
bidimetrically had failed, due largely to the tendency of these organisms to 
stratify before readings ceuld be taken. The Lumetron Colorimeter 
(Model 400-A) which is of the jewel-bearing type and has an extremely 
rapid response has proved ideal for this type of work. 


TABLE 1 
Base Mepium (Minus METHIONINE) 
1/ML. /ML. 

L-arginine 125 Biotin (free acid) , 0.0005 
L-histidine 125 Ca pantothenate 0.10 
DL-isoleucine 125 Thiamine HCl 1.00 
L-leucine 250 Nicotinamide 0.10 
L-lysine 250 Pyridoxine HCl 0.10 
DL-phenylalanine 350 Riboflavin 0.10 
DL-threonine 125 Pteroylglutamic acid 0.01 
L-tryptophane 50 Choline Cl 1.00 
DL-valine 125 Hydrolyzed yeast nucleic acid 100.00 
DL-serine 250 Factor II (from Cerophyl) 1:5 ~—_.......... 

Dextrose 1000.00 

M9S0,.7H20 100.00 

MgSO; 

K2HPO, 100.00 

CaCly.2H2O 50.00 

FeCl;.6H.O 1.25 

MnCl..4H2O 0.05 

ZnCl, 0.05 


All organisms in experimental series were grown in 4-ml. volumes in 
14.5 X 125-mm. Pryrex tubes. The tubes were incubated in a slanted 
position. The increased volume (2-ml. volumes were previously used) was 
necessary for use in the colorimeter and the slanted position increased the 
surface area for oxygenation, thereby increasing the yields. All experi- 
ments were carried through three serial transplants, readings being taken 
only on the third transplant. Our method has been to take readings on all 
third transplant tubes before inoculation, using distilled water as a blank. 
All readings are made using red filter 650. After inoculation and incu- 
bation a second reading is taken. The optical density of the un- 
inoculated tube is subtracted from the optical density obtained after 


growth has occurred. This procedure is necessary, since the media used 
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Dose response curves. Growth measured in optical density is plotted against amounts 
in y per ml. of medium. Curve 1 represents the response to 407 of homocystine per 
ml. and graded amounts of methionine. Curve 2 represents response to methionine 
only. Curve 3 represents response to methionine in the presence of 120y of L.E.L. 
Curve 4 represents response to homocystine (upper scale) in the presence of 120y of 
L.E.L. (no methionine added). All values from third serial transplants after 6 days 
incubation. 
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are colored. With this method tubes need be matched only roughly for 
size. All readings are taken with the Pyrex label as a position marker. 
The base medium is given in table 1. This is the same as Medium III as 
previously described,’ omitting methionine. Factor II is added as an 
extract of Cerophyl®. 

The liver fraction used is the same as that employed by Bennett and 
Toennies® and referred to by them as L.E.L. We are indebted to these 
authors for a generous supply of this material. In this work we were not 
concerned with the absolute potency of the active material in the liver 
fraction, accordingly it was added to the medium and the whole was steril- 
ized by autoclaving for 10 minutes at 15 Ibs. pressure. Some destruction 
of the active material may have occurred. Assays of the L.E.L. and of our 
Factor II preparation for methionine and choline were carried out micro- 
biologically. Leuconostoc mesenteriodes P-60 was used for methionine 
assay'® and a cholineless Neurospora} mutant" for the choline assay. As- 
says were made on material before and after hydrolysis to preclude failure 
of detection of bound methionine. It was found that the particular prep- 
aration of Factor II used in these experiments contained 2 of methionine 
perml. The amount used brought in 0.8y per ml. of medium. The L.E.L. 
brought in 0.4y per ml. of medium. These amounts were included in cal- 
culations when methionine additions were made. In no case was the 
amount of methionine present lower than 0.8y per ml. due to the necessity 
of adding the growth factor and where L.E.k. was also included there was 
never less than 1.27 per ml. of medium. 

Results.—It was found that homocystine could not replace methionine for 
Tetrahymena when the Factor II preparation was obtained from plant ma- 
terial. No growth resulted when methionine was withheld. This result 
was different from the one previously reported*® when Factor II was pre- 
pared from Liver fraction L. The addition of varying amounts of L.E.L., 
however, resulted in growth where methionine was replaced by homocy- 
stine in the medium. The minimal concentration of the L.E.L. for moder- 
ately good growth (under our conditions) was found to be 1207 per ml. of 
medium and this amount was used throughout. 

Dose response curves were constructed for homocystine. A typical 
curve is shown in figure 1 (curve 4). Increased growth was obtained until 
407 per ml. of homocystine were present. Higher amounts began to show 
inhibition. Asa result of these experiments 40y of pL-homocystine per ml. 
of medium were used routinely. 

In view of the fact that 120 y per ml. of L.E.L. made it possible for the 
ciliates to utilize homocystine in place of methionine, it was important to 
determine what effect the L.E.L. would have on growth in the presence of 
methionine. Accordingly two types of experiments were performed, one to 
determine the response to methionine alone and the other the response to 
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methionine plus 1207 per ml. of L.E.L. The dose response curves are 
shown in figure 1 (curve 2 for methionine and curve 3 for methionine plus 
L.E.L.). It will be noted that the minimum methionine present in curve 
2 is 0.8y per ml. (from the Factor II) while the minimum for curve 3 is 1.2 
per ml. (Factor II and L.E.L.). These curves show that the addition of 
L.E.L. with methionine makes no significant difference. 

Cystine must be synthesized, under these conditions, from methionine 
and the addition of cystine had earlier been shown to spare methionine.® 
If this synthesis takes place via homocysteine formation as suggested by 
du Vigneaud” then the addition of homocystine to the medium should 
likewise spare methionine. This was found to be the case as shown in 
figure 1 (curve 1). The addition of 407 of DL-homocystine per ml. spared 
approximately one half of the methionine. This addition made the 1.2y of 
methionine, inevitably present as a minimum when L.E.L. was used, ca- 
pable of supporting considerable growth (approximately 0.100 on the O.D. 
scale). It should be noted, however, that 0.8y of methionine was insuffi- 
cient for growth even in the presence of homocystine. - It appears, therefore, 
that some of the action of the addition of L.E.L. to the homocystine medium 
is merely the bringing in of enough methionine to raise the level above the 
growth threshold. This does not fully account for the resulting growth 
when L.E.L. and homocystine replace methionine. Under these condi- 
tions, with only 1.2 of methionine present, growth is more than double that 
obtained when L.E.L. is omitted (0.210 on O.D. scale). This difference 
must be accounted for by the conversion of homocystine to methionine. 

Discussion.—The results outlined above confirm the observations of 
Toennies and Bennett *° on the rat and point to an active agent in the liver 
fraction which in some way makes the conversion of homocystine to methi- 
onine possible. Labile methyl groups in the L.E.L. appear to be precluded 
in their experiments and in ours as the levels of methionine and choline 
could not account for enough to transmethylate the amount of homo- 
cystine used{ even if all three methyl groups of choline were available. If 
the statement of McKittrick* regarding the number of methyl groups avail- 
able from choline is considered, then the amount present is even less nearly 
adequate. The rat, unlike the chick, *! can synthesize choline” provided 
methionine is present and we have found rapid synthesis of choline by 
Tetrahymena under like conditions. 

Our results differ from those of du Vigneaud, ef al. }* and Toennies, ef 
al., ** regarding choline plus homocystine. Unlike those workers we found 
that the addition of choline made no difference in the utilization of homo- 
cystine., Choline plus homocystine did not permit growth of the ciliates in 
the absence of adequate methionine or the L.E.L. agent. This fact indi- 
cates that the L.E.L. is not supplying a methyl donor but is supplying some 
agent concerned in the transmethylation process. 
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As a working hypothesis we make the following proposals to account for 
the action of the agent in L.E.L. in transmethylation; (1) In Tetrahymena 
the usual source of methyl for the synthesis of choline is methionine and 
this reaction is made possible by what we will call ‘‘transmethylase 1”’ 
(TM1). This enzyme appears to require no exogenous co-factor. The 
reaction may be diagrammed as follows, occurring in two steps: 


Serine ———> aminoethanol'® 
TM1 
Aminoethanol + methionine ———> choline + homocysteine. 


(Il) When methionine is replaced by homocystine, methionine is synthe- 
sized by the methylation of the latter, which is the reverse of the above 
reaction. In many cases one and the same enzyme is known to catalyze 
both a reaction and its reverse. The facts here do not support such activity 
on the partof TM1. This transmethylation (homocysteine to methionine) 
requires transmethylase 2 (TM2) which must be activated by an exo- 
genous co-factor (CTM2) which the animal organism is incapable of syn- 
thesizing and therefore must obtain from its environment. This co-factor 
appears to be stored in the liver of mammals, presumably concentrated 
there from the food or from the products of microsynthesis by intestinal 
bacteria. It appears to be absent or in low concentrations in higher plants. 
Therefore, only when CT M2 (from L.E.L., in this case) is added with homo- 
cystine, is TM2 activated for the required transmethylation to methionine. 
(III) To account for the source of the methyl groups for this reaction one 
must assume the synthesis or mobilization of methyl from some unknown 
precursor (CH;R), probably first to form choline and then to form methio- 
nine via homocysteine. (IV) The transmethylation of aminoethanol by 
CH;R to form choline is presumably controlled by a trarsmethylating 
enzyme also and there is no evidence to indicate whether it is TM2 or a 
different one (TM3). It cannot be TM1 because no growth of the ciliate 
results in the presence of adequate choline and homocystine but no L.E.L. 
If it is a different enzyme from TM2 then it requires a different co-factor 
(CTMS3) and this co-factor is also found in liver. This last assumption will 
become clear from the discussion of the results on the rat. The above re- 
actions may be diagrammed as follows: 


TM2 + CTM2 (from L.E.L.), 
or 
TM3 + CTM3 (from L.E.L.), 
Step 1. CH;R > choline. 


TM2 + CTM2 (from L.E.L.), 
Step 2. Choline + homocysteine > methionine 








Homocysteine ——— cysteine (sparing of methionine). 
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In the rat !~* step 2 seems to proceed without any exogenous co-factor but 
this can be explained on the basis of bacterial synthesis, storage or both. 
Both of these sources are denied Tetrahymena so step 2 only proceeds when 
CTM2 is added. In the rat step 1 is blocked by either too low a con- 
centration of CTM2 for both steps or, if a different co-factor is required 
(CTMS3), then its amount is inadequate within the body. In any case 
Tetrahymena, needing the co-factor for step 2, has a supply provided for step 
1 by the addition of L.E.L. 

The possibility of the direct methylation of homocysteine from CH;R to 
form methionine should perhaps be considered. We know from the work 
of Simmonds, e¢ al.,"° that in the presence of added choline labeled methyl] 
from choline does appear in the methionine. This need not be the only 
route of methionine formation, especially when choline, as a vitamin, is 
being synthesized. This alternative may be diagrammed as follows: 


TM2 + CT M2 (from bacterial 
synthesis + storage or from L.E.L.) 


CH:R + homocysteine > methionine 
TM1 
Methionine + aminoethanol ——— choline (for vitamin activity). 





Recently, Horowitz” has presented evidence to show that the synthesis 
of methionine by Neurospora proceeds from cysteine and each step is 
controlled by one gene. One mutant failed to grow on homocysteine but 
did grow on methionine. It would be interesting and important to deter- 
mine whether this failure to transmethylate is due to a failure on the part 
of the mutant to synthesize the necessary co-transmethylase. This could 
be determined by the addition of the active agent from L.E.L. The fact 
that none of the Neurospora mutants could utilize homocystine and all but 
the one meritioned above could grow on homocysteine is apart from the 
transmethylation problem, and concerns only the inability of the mold to 
reduce homocystine.! 

The various parts of the above hypothesis are subject to test. Plans are 
made for the concentration and possible identification of the factor or fac- 
tors in the liver extract. Factor II preparations entirely free of methionine 
and choline can be prepared so that the possible effects of traces of these 
substances can be determined. As concentration of the active agent or 
agents proceeds tests can be made to determine whether or not more than 
one co-factor exists which will indicate the relationships between the en- 
zymes controlling steps 1 and 2. 

Summary.—1. Homocystine can replace methionine in the diet of 
Tetrahymena provided some factor or factors found in liver are present. 

2. The liver factor appears not to be a methyl donor as is indicated by 
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the failure of methionine replacement by homocystine in the presence of 


choline. 
3. The liver factor can best be described as a co-transmethylase. 


4. A discussion of the mechanism or mechanisms of transmethylation in 
animals is given. 


* Aided by a grant from the U. S. Public Health Service and a grant recommended 
by the Committee on Growth acting for the American Cancer Society. 

+ The cholineless Neurospora used in this investigation was obtained from Dr. E. L. 
Tatum of Yale University. 

¢ Betaine was absent from L.E.L. according to Bennett and Toennies. 

§ Since this manuscript went to press, we have tested strain 38706 (the methionine- 
less mutant of Neurospora crassa, referred to as Me-1 by Horowitz™). Several other 
methionineless strains of Neurospora from the Yale collection (6073A, 4195a, 6279a, 
4195A(3)) were also tested. All were found to be incapable of utilizing homocysteine, 
even in the presence of L.E.L. All responded tomethionine. This indicates that the 
genetic block occurs in the synthesis of the transmethylating enzyme and not in the 
synthesis of the cofactor. We wish to thank Dr. N. H. Horowitz for strain 38706 and 
Drs. E. L. Tatum and S. Simmonds for the other strains and for the sample of homo- 
cysteine used. 
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LINKAGE AMONG GENES CONTROLLING INHIBITION OF LYSIS 
IN A BACTERIAL VIRUS* 


By A. D. HERSHEY AND RAQUEL ROTMANT 


DEPARTMENT OF BACTERIOLOGY AND IMMUNOLOGY, WASHINGTON UNIVERSITY SCHOOL 
OF MEDICINE, St. Louis, Missouri 


Communicated by M. Demerec, December 5, 1947 


Experiments in which bacterial cells are simultaneously infected with two 
different viruses of the T2 serological group have suggested that both can 
grow in the same cell,' and that interaction between them produces new 
genetic types of virus.» * Further experiments, in which the mixed clones 
of virus coming from single bacterial cells were examined, have confirmed 
these inferences, and indicate that the new types arise only in those bacteria 
in which both viruses succeed in growing freely. These experiments are 
being continued in an effort to learn something about the mechanism of 
genetic interaction, and will not be reported in further detail at this time. 
In this paper we wish to describe a new series of mutations in the bacterial 
virus T2H,' discovered in connection with the experiments mentioned 
above, which point to the operation in this virus of a system of genetic link- 
age. 

The virus T2H, like all its relatives, exists in the form of a lysis-inhibiting 
wild type (r+) and an r (for rapidly lysing) mutant phenotype, which has 
lost the lysis-inhibiting character.1 The principal new finding to be re- 
ported here is that all the independently arising r mutants of T2H that 
we have examined are genetically different. The basic observations may 
be summarized as follows. 

A genetic difference between two 7 mutants is shown when bacteria in- 
fected with both liberate viral progeny containing a measurable proportion 
of lysis-inhibiting types, whereas bacteria infected with the same number of 
particles of either mutant alone yield only 7 viral progeny. The experi- 
mental procedure used is the one-step growth experiment of Delbriick and 
Luria, with mixed multiple infection, carried out as described earlier.’ 
For simplicity of language, we shall speak of this experiment as a cross, de- 
scribed as heterologous when new types appear among the yield of virus, 
and homologous when only the parent types are recovered. 

Certain details of the method require mention. First of all, one examines 
by it not the progeny of a genetic cross strictly speaking, but simply a 
mixed population of virus growing in the same community. The crosses re- 
ported in this paper were so arranged that the yield of virus examined came 
from about 40,000 mixedly infected bacteria. Variations in yields from in- 
dividual bacteria‘ are thus left out of account. Care was taken to ensure 
adequate and equal multiplicity of infection, with respect both to titers of 
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stocks used and to adsorption to bacteria. Fortunately, all mutants were 
found to be equally well adsorbed in sensitive tests with admixed wild type 
(table 2). The effect of unequal multiplicity on the yield of segregants was 
also tested. The largest proportion of *wild type resulted from equal 
multiplicity of infection with the two parent types, but the yield was not 
markedly affected by variations between 1 :2and2:1. Variations in total 
multiplicity above about 3 of each type per bacterium had no measurable 
effect on the result. A multiplicity of about 5 of each type per bacterium 
was habitually used. Attention was also paid to the total yield of virus in 
crosses. In parallel experiments, the yields seldom differed by more than 
15 per cent, and fell between 300 and 600 per bacterium, depending in part 
on variations in the nutrient broth. Because of these variations, difficult 
to control, important comparative results were always checked by parallel 
experiments. 

The r mutants so far examined, among which 15 have been named r1, 
r2, etc., in order of isolation, are all identical with respect to type of plaque. 
Their behavior in intercrosses places them in two clearly defined classes, 
called A and B. Crosses between any pair belonging to the two classes 
yield about 15 per cent of wild type among the total yield of virus, no sig- 
nificant differences having been established. Crosses within classes yield 
proportions of wild type characteristic for each pair, falling between 0.5 
and 8 per cent. The first eleven 7 mutants to be examined were inter- 
crossed in all combinations, yielding only one example (r1) of class A, and 
no homologous pairs. Six others were examined less completely in a search 
for a second member of class A; all proved to be class B, and at least two 
differed from all the others. The second member of class A was finally ob- 
tained from an r mutant of the distinct virus T4, by crossing with T2H.? 
This mutant is called r14. Three examples of class B were discarded with- 
out numbering. 

The results of the quantitative intercrosses, summarized in table 1, pre- 
sent three significant features. (1) The amounts of wild type appearing in 
the three possible crosses among sets of three mutants of class B (e.g., 72, 
r5 and r7) are approximately additive, but generally fall noticeably short, 
so that no linear order valid for several loci can be discerned. (2) The re- 
sults of intercrosses among sets of four mutants (e.g., r2, r4, 75 and r7) can- 
not be interpreted in terms of four independent transfer frequencies, but 
indicate some type of linkage. (3) If linkage is assumed, the classes A and 
B are probably to be regarded as independent linkage groups. However, 
the fact that rl of class A gives equal amounts of wild type with all members 
of class B is not adequate proof of this independence, for the same thing is 
seen at a lower level in the crosses between 713 and the other members of 


class B. 
The different r mutants also differ with respect to rate of spontaneous 
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TABLE 1 
PER CENT OF WILD TYPE IN VIRUS YIELDS FROM BACTERIA INFECTED WITH PAIRS OF 
r MUTANTS 

-—cLass A— - cLass B 
rl r14 r13—s 2 r4 r8 r9 13 rd r7 76 
ri 0.5 14 17 16 = (15) (13) 14 15 15 15 
r14 (13) X x xX x x x 15 xX 
113 6 6 (7) x (7) (6) 7 8 
Yeo End Bug 1.5 2.0 3.9 3.0 2.4 
7% Wt (1.5) 2.6 2.9 3.5 3.2 
78°: (0:8): st) (238). 2:8). 237 
r9 (0.9) (1.8) (2.2) (2.4) 
r3 (1.2) {E.@ 1.4 


r5 0.8 (1.1) 
r7 ~— (0.5) 
The crosses indicated by X have not been made. Results in parentheses are from 
single experiments; the others are averages of several. Self-crosses yield less than 0.1% 
wild type, usually none. 


reversion to lysis-inhibitor during serial passage of the virus.!. In two or 
more independent tests starting from different single plaques of each of 
several r mutants, the number of passages required to yield lysates contain- 
ing 10 to 20 per cent oi lysis-inhibiting virus varied as follows: for r9, two or 
three; for 75, r6, 711, five; for r2, r3, r4,°r8, 710, 14, six or seven; for 77, 
r13, eight; forrl, ten; and forrl5, about twenty-two. The lysis-inhibitor 
recovered from 72, 713, and 715 differed from authentic wild type in ap- 
pearance of plaques, being of intermediate character between r and wild 
type. On further serial passage, variants apppeared which resembled wild 
type more closely. Both types of lysis-inhibiting virus obtained from these 
three mutants were genetically different from wild type, for when either 
was crossed with wild type, excessive numbers of 7 mutants were found 
among the progeny. This was true of seven independently arising rever- 
sions from one or another of the mutants named. Presumably the lysis 
inhibiting virus selected in these cases arose by suppressor mutations, 
rather than by back mutation at the respective 7 loci. On the other hand, 
 lysis-inhibiting viruses recovered from each of the 11 remaining r mutants 
mentioned above were indistinguishable from wild type both phenotypically 
and in back cross. We have not, however, done the careful experiments 
which might reveal suppressor mutations closely linked to an r locus. 

Two different r loci were found to mutate independently in a double- 
mutant virus obtained as described in the following section of this paper. 
In four independent serial passages starting from different single plaques of 
the double mutant r1r5, the reversion passed through the more stable type 
rlineachcase. This suggests that the difference in the rate of reversion of 
the two single mutants is due to differences in stability of loci, rather than to 
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different selection rates among mutants growing in competition with wild 
type. Direct estimations of selection rates have not been made, however. 

In contrast with the spontaneous reversions, lysis-inhibiting virus arising 
in crosses between 7 mutants is always homologous with wild type. Tests 
of this point included segregants from the cross 72 X r13, two mutants 
whose spontaneous reversions are heterologous with wild type. There is, 
furthermore, no correlation between the individual stability of two r mu- 
tants and the amount of wild type they throw in crosses. 

A few of the mutants, including two lines each of r2 and 713 mentioned 
above, were carried through the transformations r — hr — r by crossing 
with a host-range mutant (1), followed by a cross with wild type. The 
resulting stocks showed the same behavior in crosses as the originals. It is 
unlikely, therefore, that these are multiple-locus mutants. 

In order to compare the ability of the different r mutants to grow com- 
petitively in the bacterial cells, mixed infections with ry mutant and wild 
type were made. In only two cases did the proportion of wild type in the 


TABLE 2 
COMPETITION BETWEEN ¢ MUTANTS AND WILD TYPE IN MIXED INFECTION 





% wip TYPE————. 





— % WILD TYPE 





MUTANT INPUT ADSORBED YIELD MUTANT INPUT ADSORBED YIELD 
rl 54 59 53 r4r7 ve Ce 55 53 
r5 50 57 54 r2 68 fief 64 
rird 53 55 37 r3 58 58 60 
113 51 48 60 16 60 60 57 
r4 60 64 65 r2r3 63 70 67 
r7 59 54 61 r2r6 60 56 64 
r4r7r13 59 54 54 r3r6 66 71 69 


yield differ appreciably from the proportion in the infecting mixture, or 
from the proportion in the fraction adsorbed to bacteria. With the mutant 
r13 the yield contained a proportion of wild type in excess of that adsorbed 
amounting to about 10 per cent of the total virus (3 experiments). The 
significance of this is not clear. The other exception is easily explained. 
The cross between the double mutant 7175 and wild type yielded an in- , 
creased proportion of r mutant amounting to about 15 per cent of the total 
virus (2 experiments), owing to the segregation of the two 7 loci described 
below. For other multiple mutants studied, all of class B, the loss of wild 
type was too small to be detected, presumably because of linkage between 
the loci. The data for individual experiments are summarized in table 2. 

Segregation of Loci.—The cross between any two r mutants yields two new 
viral types: wild type and an r mutant distinguishable from others only by 
back cross to the parental types, in which it proves homologous with both. 
This behavior supports the idea that each mutant represents a modification 
at a different genetic locus among many controlling inhibition of lysis. 
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For any pair X and Y, the segregation may be written: (X*) X (*Y)= 
(++) X (XY). In testing this hypothesis, it is necessary to prepare single- 
plaque stocks from a number of the r progeny of the cross, and test each one 
by back cross to each of the parent types. Tests of this kind have confirmed 
the hypothesis in every case. It was also found, within the limits imposed 
by the small scale of the experiments, that in each cross the two recom- 
bination types appeared in equal numbers, and the two parental types in 
the same proportion with which the bacteria were initially infected. These 
data are shown in table 3. In one experiment, virus yields from single bac- 
teria‘ infected with r4 and 77 were found to contain a small proportion of 
wild type in nearly every burst, the average being about 4 percent. From 
one of the yields, tests of a number of r progeny revealed the three ad- 
ditional viral types expected. 


TABLE 3 
SEGREGATION OF r Loci IN INTERCROSSES 
WILD TYPE e 
Cross (% OF YIELD) DISTRIBUTION OF TYPES AMONG Ff PROGENY 

rl X 75 15 4rl 475 2 rir5 

rl X17 15 9 rl 5 17 6 rir7 

rl X r13 15 5 rl 3 713 4 rlrl13 

r4 X r13 6 12 r4 17 713 1 74r13 

17 X r18 7 7 ¥1 7 713 1 r7r13 

12 X 73 3.3 36 12 47 r3 3 r2r3 

12 X 76 3.4 33 12 41 r6 1 r2r6 

r3 X 16 1.9 63 713 23 16 1 r3r6 

++ X rlri 37 2 rl 2 75 6 rlrd 

rir5 X rirl3 0.0 13 r1rd 13 r1r13 2 rlr5r13 271 
r4r7 X r7r¥3 0.0. 12 r4r7 15 r7r13 3 r4r7r13 017 


The equality of numbers of the two recombination progeny, only sug- 
gested by the data of table 3, has been confirmed on a statistically adequate 
scale for crosses between host-range and r mutants to be reported elsewhere. 

Multiple Factor Tests —In an effort to learn something about the nature 
of the linkage, multiple factor crosses have been made. These are of three 
types. 

The first type consists of three factor crosses involving the loci rl, 74, and 
r7, the two latter being rather closely linked but probably independent of the 
first. If this is correct, one would expect the factor 71 to distribute itself 
equally between the two types of virus resulting from the enchange of 74 
with respect to 77. The effect should be an equal segregation of wild type 
from the crosses rlr4 X r7 and rlr7 X r4, amounting in either case to about 
half that segregating in the cross r4 X r7. For the cross rl X r4r7, the 
yield of wild type should be nearly as large as that for rl X 77, owing to the 
linkage. The data shown in table 4 confirm these expectations very well. 
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TABLE 4 
SEGREGATION OF WILD TYPE IN MULTIPLE FACTOR INTERCROSSES 


CROSS % WILD TYPE CROSS % WILD TYPE 
rir4 X 17 1.4 (8) r7r13 X r4r13 X r4r7 0.16 (2) 
rir7 X r4 1.4 (3) r4 X r7 X 713 7.0 1) 
rl X rAr7 11 (2) r2r3 X r2r6 X r3r6 0.02 (7 
rl X17 15 (4) 12 X 13 X 16 3.0 (2) 
r4r713 X 17 1.8 (4) 727r3 X 16 0.8 (7) 
r7r13 X r4 0.7 (3) r2r6 X 73 0.3 (6) 
r4r7 X 113 5.0 (3) r3r6 X r2 1.4 (5) 
r4 X 713 6.3 (5) 12 X 73 2.0 (6) 
v7 X 713 7.0 (8) 12 X r6 2.4 (6) 
r4 X 17 3.5 (9) r3 X 76 1.4 (7) 


The results shown are averages of the number of experiments indicated in parentheses. 
The 22 measurements, including all crosses r4 X 113, r7 X 113 and r4 X r7 made during 
several months, show a coefficient of variation of 20%. 


The second type of cross corresponds to the conventional three factor test 
for linkage order. We have used two sets of three factors all belonging to 
the class B. The sets include (72, r3 and r6) and (r13, r4 and r7), the 
former being more closely linked than the latter. The data for these 
crosses, also included in table 4, tend to place the factors in the order listed 
above, though the correspondence with expectations for random crossing 
over between linear structures is rather poor. The data for r2, 73 and r6, 
particularly, are compatible with the possibility that a crossover in one 
region tends to be accompanied by a second in the adjacent region. 

The third type of cross, suggested by analogous experiments with bac- 
teria,® was tried for the purpose of estimating the amount of repeated ex- 
change (i.e., exchange between recombination progeny and progeny of the 
parental types) taking place under the conditions of the experiments. The 
combinations used were such that segregation of wild type would require 
the interaction (simultaneous or successive) of three different mutants. 
The sets chosen were (74713, r7r13,r4r7) and (r2r3, r2r6, r3r6). In dif- 
ferent experiments, multiplicities of about three and five per bacterium of 
each of the three mutants of the set were tested, with similar results. The 
data for these crosses, and for the corresponding sets of three single mutants, 
are shown in table 4. It will be seen that repeated exchanges can be de- 
tected by this method; that they are less frequent for the more closely 
linked factors; and that their number is too small to be of quantitative im- 
portance in the linkage tests already described. 

It should be noted that repéated exchange of the type detected above 
should not have any effect in two factor crosses. Repeated exchange be- 
tween sister recombination progeny, if such arise, presents a more serious 
difficulty. Until effects from this cause can be evaluated, no great quanti- 
tative significance can be attached to the results of intercrosses. 
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Discussion.—Differences among the r mutants described in this paper 
indicate that there are an indefinite number of distinct loci controlling lysis 
inhibition in the bacterial virus T2H. The suggestion from the genetic 
data is that the mechanism of lysis-inhibition®’ may be exceedingly complex. 
On the other hand, the fact that identical mutants are not encountered 
shows that the different ry mutations occur with similar frequency. This 
remarkable circumstance might be interpreted as an indication that the 
different mutations result from a single type of chemical change occurring 
at different points in the genetic material of the viral particle. If this were 
so, the different ry mutants might represent alterations of the same biochemi- 
cal property. This argument must be qualified in view of our ignorance of 
the number of possible 7 mutations. 

The existence of suppressor mutations noted in this paper, and of inter- 
mediate mutant types,* complicates the system somewhat further. 

It is evident that the rather high frequency of r mutations, about 10~‘ per 
viral duplication,’ is at least partly due to the large number of loci involved. 
It is probably not feasible to measure the stability of any given r locus. 
For this reason the test of equal] rates of y mutation in wild type and in the 
identical phenotype obtained by reversion from 7,* is of no value as a test of 
reversibility of mutation. The test by crossing with wild type, reported in 
this paper, shows that some reversions are not back mutations, but that 
most of them probably are. 

The situation described for the 7 mutants does not prevail among viral 
mutants in genéral. We have examined six independently arising host- 
range mutants of the h* phenotype,* which proved to be homologous in in- 
tercrosses. This locus can be assigned a position in the linkage class B how- 
ever, and it is for their usefulness in mapping that the y mutants are likely 
to prove most valuable. . 

The linkage system itself we shall discuss only briefly, pending the com- 
pletion of experiments of another type. It is clear that the clases A and B 
are independent, or nearly so, and that some kind of linkage exists within 
the classes. Regarding the linkage structure, some of our data suggest 
linear arrangement with crossing over, but others are somewhat contradic- 
tory. Among difficulties which cannot be discussed very intelligently at 
present, the necessity for reconciling our data with the results of Luria’ 
should be mentioned. It is probably necessary to add that the linkage 
system discerned in the virus may yet prove to have little in common with 
the more familiar linkage systems in cellular organisms. 

Summary.—Different r mutants of the virus T2H differ (a) in rate of 
spontaneous reversion to lysis-inhibiting types; (0) in the kind of reversion 
they undergo, either to a form genetically like, or to forms unlike, the wild 
type; (c) in the amount of wild type they throw in intercrosses. With re- 
spect to (c) there are two independent linkage groups of distinct r loci. 
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Analysis of one of these groups by single and multiple factor intercrosses 
has failed to establish or refute the idea of crossing over between linear 
structures. The existence of some kind of linkage system conditioning 
segregation and reassortment of genetic factors among viral particles seems 
quite clear, but it remains to be learned whether reciprocal, or even ma- 
terial, exchanges are involved. 


* Aided by a grant from the United States Public Health Service. 

{ On leave from the Instituto Bacteriologico de Chile, Santiago, Chile. 
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A THEORY OF MINIMAX* 


By Max SHIFFMAN 
NEw YorK UNIVERSITY + 


Communicated by S. Lefschetz, January 7, 1948 


1. Variational analysis, or the theory of critical points, can be examined 
either analytically or topologically. To furnish the link between these 
two aspects, it is necessary to establish certain basic properties of the 
functional in question beyond what is called for in the case of the pure 
minimum. Lower semi-continuity and compactness are adequate for 
minima. But precisely this semi-continuity, and not continuity, introduces 
new questions in the more general case of minimax. A link between the 
analytic and the topologic is established if the functional can be shown 
to be reducible. This link is completed if it can be shown that the minimax 
type defined topologically is equal to the minimax type defined analytically. 
Only when these facts are established is it possible to pass between the 
domains of analysis and of topology. 

In the present note, we shall sketch a theory which does not separate 
these two basic questions, reducibility and type, but rather makes an 
organic unity of them (Theorem 3). We discuss a single integral problem 
where the functional ranges over curves, but it will be clear that the theory 
is in principle applicable to problems involving surfaces or higher dimen- 
sional varieties. Indeed these latter problems are the main motive and 
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goal of the present work and will be considered-in the detailed publication. 
2. Consider the variational problem for 


Ily] = Sea’ fo, 9, 9’ )dee (1) 


with prescribed end-points y(x0) = yo, y(%1) = yi. Let y = o(x) be an 
extremal passing through these points. We shall suppose that f is suffi- 
ciently often differentiable in its arguments x, y, y’ and that 


Fury (x, y, p) >0 (2) 


for all y in a neighborhood of ¢(x) and for all p, i.e., we suppose that (1) is 
a regular variational problem near y = g(x). The second variation of 
I{y] is 


82 [n] = SS [fan + fen’ + fren’ \dx (3) 
where the arguments in the terms involving f are (x, o(x), y’(x)) and where 
n(xo) = n(x) = 0. (4) 


The quadratic functional (3) subject to any normalizing condition, say 
JS n’dx = 1, has certain successive eigenvalues and the corresponding 


eigenfunctions. Select the first of them, di, Az, ..., A, and (x), u(x), 
.., U,(x), respectively, where 1 is chosen so that 
An+1 > 0. (5) 


The eigenfunctions u;(x) satisfy 
d 
(fyi + frm’) ine ax € fyrylt + fyryui’) hors Nii = () 


U(X) = Ui(x1) = 0 (6) 
and are normalized so that 
Sa udx)u;(x)dx = by. (7) 


With these preliminaries, the following fundamental theorems of the 
theory can be stated (for Theorems 1, 2 a much weaker hypothesis than (2) 
can be made). 

THEOREM 1. There is a family of functions 


y = o(x; C1, ...5 Cn) (8) 


depending on the n parameters , ..., Cn and defined for sufficiently small 
absolute values of these parameters, which lie in some neighborhood in the 
(x, y) plane of y = (x), which reduces toy = g(x) fore, = ... = Ge = 0, 
and which passes through the end-points {xo, yo), (x1, yi). For any curve 
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7 = 


y(x) passing through’ the end-points and lying in this neighborhood of 
y = g(x), determine the values of c, ..., Cn by 


Se (yx) — oe) u(x)dxe =c, i=1,...,0 (9) 


Then 
I[y(x)] 2 I[e(x; c)]. (10) 
THEOREM 2. For the family o(x; c) we have 
I[e(% ¢)]; = I[e(e)] + Lifer dc? + 7 (11) 


where the remainder term r is estimated by 


[r| 7 
—ASitns: sg tis 0 
oe +... t+,” : 
uniformly in C1, ..., Cn 
THEOREM 3. Make the additional assumption that there are positive 
numbers A, M such that 





| fy(x, y; p) | < 
——————| S$ M for all =A (12) 
Pfy (x, y; p)| f lel 
and for all (x, y) in a neighborhood of y = o(x). Consider the quantity 
I(t) = I[p(e; c) + tly) — ox; o)}] (13) 
as a function of tin 0 < t < 1, where cq, ..., Cx are determined from (9). 
Then 
dI(t) a“ 11) — I(0) aye 
ee t ¥ 
and 


I(t) S$ 10) + #{I(1) — I(0)}. 


3. To prove Theorems 1, 2, consider the problem of minimizing I[y(x) ] 
among curves passing through (xo, Yo), (*1, 1) and satisfying the side condi- 
tions (9) for fixed but small values of c, ..., Cn. The extremals of this new 
variational problem are solutions of 


d n 
fy — Po Sun yore 1 mates(x) = 0 (14) 


where i, ..., fn are constant Lagrange multipliers, and they must also 
satisfy (9). Let 


V(x; b, Biy e+ey Hn) (15) 
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be solutions of (14) with initial conditions 


Y (xo, b, bh) ==. Ve 
V(x, b, uw) = b+ Dofer 5 ui’ (x0). (16) 


A particular solution of (14), (16) with m, = ... = u, = Oand b = 9! (x0) 
is o(x). 
We establish first that 


Y, (x, 9’ (2%), 0) = Ce). (17) 


Then insert (15) for y(x) in (9) and consider the resulting equations as ” 
equations for the determination of mw, ..., u, in terms of b, Gq, ..., Cn. 
For b = ¢'(%), 4 = ... ='Cn = 0, a solution of these equations is 4. = 

. = pn = 0. By virtue of (17), we prove that the Jacobean of these 
equations, evaluated for b = ¢’(%), 1. = ... =, = 0,is #0. Therefore 
we can solve the equations for mw, ..., un in terms of Db, , ..., Cn for b 
sufficiently near y’(xo) and G, ..., ¢C, sufficiently small in absolute value. 
Substituting in (15) gives a function 


Z(x, b, Ci, «oy Cn) (18) 


which is an extremal and satisfies (9). 
We establish now that the equation 


Z(x1, b,c) = 1 


has a unique solution for d in terms of ¢, ..., Cn for sufficiently small values 
of these parameters. Substituting in (18) gives the final function 


y = o(x; G1, ..., Cn) (19) 


which is an extremal, passes through the given points and satisfies (9). 
For the extremal (19), the second variation of the new variational 
problem is 


- { Fn? + fyynn’ + Fyn’ } dx (20) 
with 
Sun nu(x)dx=0, t=1,..., 0, (21) 
and (4). The bar over the quantities involving f in (20) indicate that the 
argument is (x, 9(x; c), y’(x; c)). The quadratic functional (3) subject 
to (21) and (4) is positive definite, and has a value 2 An+1 JS n'dx. Since 
v(x; c) and g’(x; c) are near y(x), (x), the second variation (20) subject 
to (21) and (4) is also positive definite. The positive-definiteness of the 
second variation implies that the extremal (19) can be imbedded in a 
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suitable field of extremals. Furthermore, by virtue of the condition (2) 
and the fact that the side conditions (9) do not involve y’, the Weierstrass 
E-function for this variational problem is positive. Therefore (19) 
actually minimizes, and this proves Theorem 1. 

Theorem 2 follows from the method of determination of g(x; c). For 
we prove that ° 


Po,(X; 0, oo 8p 0) = u(x), 


and this will yield Theorem 2. 

Incidentally in (16) we tacitly assumed that none of the \; are zero. If 
some \; = 0, a suitable modification of (16) can be made. 

4. To prove Theorem 3, consider the integral 


J{w] = Ens g(x, w, w’)dx (22) 
where ; 


g(x, wy w’) = wf, (x, ¢(x; c) + W, g’ (x; c) + w’) 
+w'fy (x, o(x; 6) + w, 9 (x; 6) + o’) (23) 
—Yef(x, o(x;c) + w, 9'(x;c) + ow’) 


and where w(x) is subject to the side conditions 
Sn wul(x)dx =0, i= oe n (24) 
and the end conditions 
; (x0) = w(x) = 0. (25) 


The variational problem for J[w] subject to (24), (25), which we shall call 
problem J, has as Euler equation 
d 
ee. > eee 2 
a7) dx Su ane Hits(2) 0. ( 6) 


But for the curve w = 0, we have 


fo = V/efy, be = Vofy 


where the bar over the f indicates that the argument is (x, y(x; c), 9’ (x; 


c)). Since g(x; c) satisfies (14) for suitable values of i, ..., un, it follows 
that w(x) = 0 satisfies (26) for these constants mw, ..., Hn and so 
w(x) =0 (27) 


is an extremal for problem J. 
Evaluating for the extremal (27), we have 


80'e’ = 8/otyy > 0. (28) 


The second variation of J is likewise */, times the expression (20), with 
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side conditions (21). The extremal can therefore be imbedded in a field 
of extremals for the problem J. 
The Weierstrass E-function 


E= g(x, , g) ssn g(x, ®, w’) igs (o kt wo’) gur(x, Q, w’) (29) 


for w, w’ in the field, and o arbitrary, can be shown to be positive in a 
sufficiently small neighborhood of (27). For from (28), there is a constant 
a such that E > 0 for |o| < a(ando #w’). By virtue of (2) and (12), 
there is a constant B such that E >0 for |o| => B. Finally, in the range 
a < |o| S B we can make E > 0. In all three cases, sufficiently small 
neighborhoods of (27) in the field must be taken. This establishes the fact 
that (27) actually minimizes J[w] in the final neighborhood selected. 
Hence 


J[w(x)] 2 —'/el[e(x; c)). (30) 
Consider the quantity J(t) defined in (13). We have 
dit) lay _ 
—_ ~ 5 1) = J[w] (31) 


where 
w = tHy(x) — o(x; 0)}. 
Inequality (30) yields 


dit) _ 1 1 

in _- 9 T(t) 2- 9 1(0). 

This is the first inequality of Theorem 3. The second inequality is an 
immediate consequence of the first. 

5. Theorem 3 can be formulated as a fundamental deformation 
theorem which establishes the reducibility of J[y] and the characterization 
of the minimax type of the extremal y = g(x). Select m so that A,+1 is 
the first positive eigenvalue. There is a neighborhood of y = g(x) such 
that any curve y(x) in this neighborhood can be deformed into ¢(x; c) 
in such a way that the functional J[y] always decreases, and at a rate 
given in a uniform way in Theorem 3.. The end family ¢(x; c) is an n- 
parametric family for which Theorem 2 holds. This proves reducibility 
and at the same time shows, in case none of the ),; are zero, that the topo- 
logic type of minimax is exactly equal to the number of negative eigen- 
values of the associated boundary value problem (6). 


* Dedicated to Richard Courant on his 60th birthday. 
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ON ALL OF MERSENNE’S NUMBERS PARTICULARLY Myo 
By H. S. UHLER 
YALE UNIVERSITY 


Communicated by F. D. Murnaghan, January 15, 1948 


On November 27, 1947, the author finished investigating, by application 
of the Lucasian sequence 4, 14, 194, 37634, ---, the factorizability of 
Mersenne’s number Mj, = 2% — 1 = 1255 42034 70773 36152 76715 
78846 41533 28322 04710 88892 80690 25791. The 192nd residue had the 
value 542 45701 25193 90814 13211 43009 56802 04633 04970 79432 42801 
51282 which, being non-zero, shows! that Mj; is composite in-character. 
A sufficiently detailed account of the procedure followed and of the under- 
lying basic theorem may be found in the author’s first paper? on Mersenne’s 
numbers. Suffice it to record two facts. After the date given above the 
. data of the entire set of original work-strips were checked for the third 
time with an auxiliary modulus in conformity with the condition (Migsqx + 
r; + 2) = r—2(mod 10’ + 1), where g and r denote “quotient” and 
“remainder,’’ respectively. The approximation to the reciprocal. of 
M3 used in the latest work had the value 0. (58 zeros) 79654 59555 66226 
13851 44401 98883 85590 27955 52277 59630 93930 37006 37523 90143 
41987 11093 10854 40700 94876 19334 41734 25085 25958 34273 96290 
7--++. Thecheck multiplication gave (Mi3)(1/Miss)a = 1+ 1.98 K 10718, 

With regard to the 55 numbers of the form 2? — 1 (where # is a prime 
not exceeding 257) considered by Mersenne, the following brief report on 
recent progress made in their study and the present state of our knowledge 
may be timely because of its definitive character. In the year 1935 there 
remained*® six Mersenne numbers which had not been investigated with 
respect to their prime or composite properties. These corresponded to 
p = 157, 167, 193, 199, 227 and 229. 

In leisure hours the present writer began the investigation of Mys7 in 
the spring of the year 1944 and finished the entire set with 4/1; on the 
date given above. These M,’s were taken up in random order in the 
vain hope of discovering a prime number greater than 2! — 1. More 
specifically the final residues for Mis7, M1s:, M229, Mig and Moo; were ob- 
tained, respectively, on Aug. 11, 1944, Dec. 2, 1944, Feb. 9, 1946, July 27, 
1946, and June 4, 1947. (The manuscript on M27 was accepted for pub- 
lication in the Bulletin of the American Mathematical Society on July 17, 
1947.) Since all six M/,'s were found to be composite the earlier status of 
Mersenne’s remarks has not been changed but a lacuna in our knowledge 
has been filled. The next extremely difficult step will consist in the dis- 
covery of all of the as yet unknown factors of Mersenne’s numbers. 

Mersenne said* that the only values of p not greater than 257 which 
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make M, prime are 2, 3, 5, 7, 13, 17, 19, 31, 67, 127 and 257. Comparison 
of this list with the correct data recorded in the top line of the table pre- 
sented below shows that Mersenne made five mistakes. p = 67 and 257 
do not yield prime values for M/,, and p = 61, 89 and 107 were not included 
in his list of special primes. 

. With reference to explicit factoring, attention should be called to a 
valuable paper’ by Professor D. H. Lehmer entitled ‘‘On the Factors of 
2" + 1.” His investigations on 76 numbers unveiled eleven factors 
which fall within Mersenne’s range. Incidentally two of his new factors 
confirmed the present writer’s final residues for Miz and Moy. 


p CHARACTER OF Mp 
2, 3, 5, 7, 18, 17, 19, 31, 61, 89, 107, 127 Prime 
11, 23, 29, 37, 41, 43, 47, 53, 59, 67, 71, 73, 79, 113 Composite and fully factored 
151, 163, 173, 179, 181, 223, 233, 239, 251 Two or more prime factors found 
83, 97, 131, 167, 191, 197, 211, 229 Only one prime factor known 


101, 103, 109, 137, 139, 149, 157, 1938, 199, 227, 241, Composite but no factor known 
257 


1 Lehmer, D. H., Jour. London Math. Soc., 9-10, 162-165 (1934-1935). 
2 Uhler, H. S., Proc. Nat. Acad. Sci., 30, 314-316 (1944). 

3 Archibald, R. C., Scripta Mathematica, 3, 112-119 (1935). 
4Lehmer, D. H., Bull. Amer. Math. Soc., 53, 164-167 (1947). 


NEW TYPES QF CONGRUENCES -INVOLVING BERNOULLI 
NUMBERS AND FERMAT’S QUOTIENT 


By H. S. VANDIVER 
DEPARTMENT OF APPLIED MATHEMATICS, UNIVERSITY OF TEXAS 
Communicated January 17, 1948 
In another paper! the writer gave the relation * 
r 
r _1 Sn(m, Rk, a 
(mb +B = (7) (-ye BO: p> as) 
a=1 


where we define the Bernoulli numbers b, by means of the recursion formula 
(b + 1)” = d,; n > 1, the left-hand member being expanded in full and 3, 
substituted for b*, and the left-hand member of (1) is interpreted in the 
same way; 


a-1 
S,(m, k, a) = xX (im+k)"; QO-=1, 


mand k are any integers withm #0. Inthe present paper we shall employ 
(1) as well as other known relations to obtain various congruences. The 
principal results proved seem to be (5), (6b), (7) and (20). 
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By the well-known Bernoulli summation formula we also have 
(b + k)" = mSy-i(k) + dn, (2) 


where 
Sn—1(R) = Reamer 0, k). 


Set in (1),7 = p — 1 we have, ifn < p — 1, p being an odd prime, 


p-1 a 
(mb + k)” = (i (? . ') (—1)2- Salm, k, a) (3) 


a=1 a 


It is anaien that 
wy ve 
(? > 1) =(-1%mod 2), (3a) 


so that from (3) we obtain 


pot 
7 Salm, k, a) 
a=1 a 


(mb + k)" = (mod p). (4) 


For a particular a = 7, the corresponding term on the right may be written 


—ljzn n ; n 
se + (k +m) +... (& + @ — 1m) ae 


t=1 4 





Now collect the terms in the right-hand member of (4) involving each 
particular nth power, we find 





p—-2 1 1 1 
C= Ee + ms" (5 + ygt tym). 





re ae p—1 
But also 
1 1 
L+ gt... +55 Hlth 7% 4+... + (6-1? % =0 (mod 9), 
(4a) 
since pis odd. The last result gives 
- (++ i +t) mitit 2 
44 s+2 sb regaig” = 2 eee _ 
so that (4) may be written, n < » — 1, 
p-2 
(mb + k)” = 2, (ma + k)"Ra(mod p), (5) 


where 
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For k = 0, m = 1, this becomes 
p-2 : 
bn = > a"R,(mod Pp). (5a) 
. a=1 
Also, (5) gives for g any integer, 0 < g < p, for mg + k # 0 (mod p), 


p-2 


DX (mg + k)~"*?-?(mb + ky” 


n=0 
ore 
= Li Re 2 (ma + k)"(mg + k)?-?-" 


(ma + k)?~* — (mg +k)?" 
m(a — g) 





p-1 
= (me + bP? Re + O Re 


axg 


or 


p-—2 
Dd (mg + k)~"t?-?(mb + k)” = —(mg + k)?-?R, (5b) 


n=0 


modulo p, by Fermat’s theorem, provided mg + k # 0 (mod p). For 


m = 1, k = 0 this reduces to 


De" = —Rymod ), (6) 
which for g = 1 becomes? 
bo + bi + ... + bp-2 = —1 (mod fp). {6a) 
Taking the relation (5a) for nm = 0,1, ..., p — 2 in turn introducing 


powers of the integers m, prime to #, similarly to the procedure in the proof 
of (5b) we obtain 


p—2 
2 e"m"b, = —Re(mod P), (6b) 


where am = g (mod p); a< p. 
Now take (5a) again and let p = 1 + cmand putn,n+c,...,." + 
(m — l1)e, in turn in this congruence, and add, we obtain if n < ¢ 


m—1 
Da batt =m >-r"R,(mod 9), (7) 
where r ranges over the distinct solutions in the set 1, 2, ..., p — 1 of 


x° =1(mod p). The relation (7) has a type of analogy with the relation® 


pag m—-1 bn 
aa a ar ~ (mod #), (8) 


k=0 
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with cand even. Set ux =k/m (mod p);0< ue< p. Then 
k n 
(mb + k)” = m” (2 ‘) =m"(b + px)” 


= m"(b, + nO" + 1"-1 + 2.2. + (ue — 1)"7")) (mod 9), 
on applying (2). This gives 


‘dh 18) Mla. an "(07 * + 18) +. + (ee — 1)*7”), 





n 
modulo p. Setting m = 1, 2, ..., p — 1, in turn we obtain by adding, 
since 

9 p—2 (ma)? * ie 1 
1 + (ma) + (ma)? + ... + (ma) = ne, 
unless ma — 1 =0 (mod ), 
p-1 — % 
SP ninth aka (9) 





i= 1 4 


according as m,; 2 mx OF w; < wy. In particular the right-hand member is 


m if k = 1. Using the known relations, 
Sn(p) = pb»(mod p), 


0 < g < p, the first congruence following involving p terms on the left, we 
have - 


gr titi+... +1) = phy’ *(mod p), 
(1 +24... +p — 1) = phig? (mod p%), 


ty ce 1)P* gt, lnod 9%), 


and we obtain by addition 


eb oe—) get R 
pe 4 Sm Beem, 
meee : 


employing (6), and if we write, with (a, p) = 1, 


a?-* — ] Z 
= 
p 
then the above gives 
p-1 ae 
ft 4 SMO) = A) _Rwog 1) 
a=1 @ — 2 g 


axe 
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Now modulo p, 


o24'5 g(a) "5 go q(g) _%~ ae ~ a(g)g?-? + gh? 





= a= ¢ — ~* ee 

using 

p—-1 —2 —2 

Sag BE — wht me mod 2), 

e=14 g a=1 

axg 
hence 

1 qe) x a) Q. —Feamoa p), (9a) 
le see ite « 


which for g = 1 gives 
°S g(a) 


Ph sper —2(mod p). 


We now obtain another formula related to (5b). Employing! (relation 


(4)) 
(b(m, k) + m)"** — dazi(m, k) = m(n + 1)k?-?, 


and setting n = p — 2, we have after employing (3a), 
p- 
> (—1)'b,(m, k)m?—1—-* = —mk?-?, 
i=0 

modulo p, where (mb + k)° = 6,(m, k). 


We shall now modify some of the ideas in a certain proof* of (8) to 
obtain new types of congruences involving the Bernoulli numbers and 
Fermat’s quotient. We have‘ 


1 


(1 - moe os 


nm” 


kp/m 
pe "(mod P), (10) 


where [x] is the greatest integer in x, m >1,n > 1, (m, p) = 1. For’ 
= 1 we shall now show that 


m—1 [kp/r] 


(m — 1), =>> x i"—'(mod p). (10a) 


fe have for (k, m) = 1, (m, p) = 1. 


ph = ret m{ O<r<m. 
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Sethi 2 3); . ae ~ 1 in turn and add we have since pk, = pko(mod m) 
only if ki = k: and therefore the r’s range in some order over 1, 2, ..., 
m — l, 


m—-1— 
as | a 7 ANS = 0 (mod p), 
k=1L™ 





a 2 
or 
m—1f- 
x eh | = b,(m — 1) (mod p). 
k=1L™ 
Now take (10) with (10a), put m = 1, 2, ..., p — 1, in turn and add, after 
multiplying each relation in turn by g'1—", we find 





p-1 p-2 s 

m —1 (m — 1) -1-s;l—m 

eee? paeel T Aees +z" So, me 
p P pra 


m 2 s=2 
—g’~? a(g, m)(mod p), (11) 


since in (10) and (10a) when we add we have on the right 


m—1 [pk/m] = ; *» ‘ ee 
DL ee" + ig? + eh +... + HP), 
k=1 t=1 
or, if 0 < g< fp, 
m1 [kp/m] gp-1 _ gh—1 a 
+” — Dele, m), 
k=1 i=1 aes sy 4 
tg 
where g®~? occurs a(g, m) times in the right-hand member of (10) and 


(10a). Now (11) is satisfied identically for m = 1. It may also be 
written 


a mol P-2,,b-1-s _ 4 eet aa 
(m 1)bp-1 : + pa : g b, = 
* poepy 
se ported m)(mod p) (11a) 
Set g = 1,m = 1, 2, ..., p — 1 and add, we have, since for m < p, that 
a(1, m) = m — 1, 
5. 
(Sps(b) — B+ bor + D0 (mod p), (12) 


since 


Sa(p) = 0 (mod 9) 


fora < p — 1. Now we can write 


Sp-i(b) = pbp-s(mod 9’). (13) 
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Using (13) and pbp_1 = —1(mod ), together with® 


— 1)! 
bp-1 =.1 + P= (nod p), 


we obtain the known® formula 
D, & = Wimod 0), (14) 


where JV is the Wilson quotient 


O31 > 7. 
p 


On using this in connection with (lla) for m < p we have another known‘ 
result 


p—2 b, 
—q(m) + 2) —, = W(mod p). (15) 
s=1 5M 


Multiplying this last congruence through by a‘ gives,® if ¢ is in the range 
|i dar o. »~P — 2, 


o-1 
YD a'g(a) = — 70d 2). (16) 


The relation (15) is an analogous theorem to (6). If m > 1, (m,p) = 1, 
p prime and odd, then it is known that’ 


late 1 Dn n— _ 
(m" — bn =>’ = yet pa (mo d p) (17) 
where 
fale) = OF + the + Me? + 20. + (p — 1)*x?? 
with 0° = 1. In this relation set n = 2, 3, ..., p — 1, together with the 
relation 


oe ptpert... +! 
phe 3 11>, 





and add these p — 1 congruences, we obtain, if r < p, the symbol 2 thy 
indicating summation over all the distinct mth roots of unity excepting 
unity, 


See + Fst... + 5?-*) 





s=1 : 1 — p? 





SR DER ni RE Das Se NBO IE ah Tt 


abe 
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r?—2(5 — Ip" et gt es 
=P DT - ot) 
Ae Soak a, a 
ye? 
=l5-4 et 


modulo ~, so we obtain, if we set 





B(m, i) = Lo. 


— 1) Z 
form >1,0<r< bp, 
=? 5 aa m" p-l—n =? 
—q(m) + 2d ar bn =r ~ B(m, r)(mod p). (18) 
But by (15) we also. have 
Ce 2 -—2 fe 
——— -—gr)+ > = W(mod p). (19) 
Now (18) may be written 
p-2 fe p-2 mr? —1—"p, 


—q(m) + 2» = r°~*B(m, r), 


n=2 n 


modulo p, which with (19) gives, ifm >1,0<r< pp, 


p—2 p— n p—1—n 
—am) + a) + + W- oe =r?*B(m, 7), (20) 


modulo p. Forr = lin Samay this takes the form 


2 m” 


= B(m, 1). (21) 





n=2 


Note that there is a certain type of analogy between (6b) and (20) with 
b, corresponding to b,/n. 


1 Duke Math. Jour., 8, 576 (1941). ; 

2 This relation is certainly one of the simplest the author has ever noted concerning the 
Bernoulli numbers. I have failed to find it in the literature, but do not imagine it 
could be new. Another proof may be obtained of it by the use of (2) by setting » = 
pb — 1,and k = p — 1 therein. 

3 These PROCEEDINGS, 31, 59, 312-314 (1945). 

4 Annals Math., 18, 114 (1917); Duke Math. Jour., 3, 571, relation (7) (1937). 

5 This relation is well known; for a recent proof, cf., Vandiver, Duke Math. Jour., 8, 
578 (1941). 

6 Lerch, Math. Annalen, 60, 471-490 (1905); Nielsen, Oversigt Danske Vidensk. Selsk. 
Forhandlinger, 518-519, 177-180 (1915). 

7 Included in a result due to Frobenius, Berlin Sitsungsberichte, 827 (1910). The 
latter theorem was generalized by Vandiver, Trans. Amer. Math. Soc., 51, 515 (1942). 
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(RITICAL POINTS OF HARMONIC FUNCTIONS AS POSITIONS 
OF EQUILIBRIUM IN A FIELD OF FORCE 


By J. L. WALSH 
DEPARTMENT OF MATHEMATICS, HARVARD UNIVERSITY 


Communicated January 29, 1948 


The first mechanical interpretation of the critical points (zeros of the 
derivative) of a polynomial p(z) is due to Gauss,! who showed that the 
critical points not multiple zeros of p(z) are the positions of equilibrium in 
the field of force due to unit particles situated at the zeros of p(z), where 
each particle repels with a force equal to the mass divided by the distance. 
Bécher (1904) showed that the zeros of the jacobian of two binary forms 
and of certain other algebraic invariants can be expressed as the positions 
of equilibrium in a similar field of force due to positive and negative particles 
situated at the zeros of the ground forms. Bdédcher’s results can be applied 
in the study of the critical points of rational functions. 

A given harinonic function can frequently be expressed as the uniform 
limit of the logarithm of the modulus of a variable polynomial or more 
general rational function; the critical points of the harmonic function are 
the limits of those of the rational functions; the present writer (1933 and 
later) has thus applied the extensive theory of the location of the critical 
points of a rational function to the study of the critical points of a harmonic 
function. It is the object of the present note to show that the latter critical 
points can be conveniently and effectively studied directly as positions 
of equilibtium in a suitably chosen field of force, due to matter spread over 
Jordan curves or other point sets. The advantages of the present method 
are a closer analogy between the theory for harmonic functions and that 
for rational functions, with simpler proofs and stronger theorems, and 
in addition a systematic method for the enumeration of critical points of 
harmonic functions. We include also the means for studying various new 
harmonic functions. 

A finite critical point of an analytic function f(z) is a zero of the derivative 
f'(z); the order of the critical point is the order of the zero of f’(z). The 
point at infinity is a critical point of an analytic function f(z), and of order 
k, if z = Oisa critical point of f(1/z), and of order k. If f(z) has a pole at 
infinity, that point is not a critical point; if f(z) is analytic there, the point 
at infinity is a critical point of order & if and only if f’(z) has a zero there 
of orderk + 2. A point z = x + ty isa critical point of a harmonic func- 
tion u(x, y), and of order 8, if it is a critical point of the analytic function 
u + iv, where v is conjugate to u, and of order k. 

We shall consider critical points as positions of equilibrium in a field 
where the force is the conjugate of the derivative f’(z) of an analytic func- 
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tion f(z). Corresponding to the convention just introduced, we consider 
the point at infinity to be a position of equilibrium if and only if it is a 
critical point of f(z). 

1. We use the term Jordan configuration (of the extended plane) to 
denote either a continuum composed of a finite number of Jordan arcs or 
a finite number of such continua, mutually disjoint. Various integrals of 
potential theory representing a harmonic function u(x, y) in a region R 
may be taken over the boundary B of Rif Bisa finite Jordan configuration; 
ordinarily such an integral is first taken over a level locus Z with no multiple 
points, say u(x, y) = const interior to R, and on L involves (Ou/dv)ds = 
—dv, where v(x, y) is conjugate (locally single valued) to u(x, y) in R, 
and v represents interior normal. If u(x, y) is constant on a Jordan 
configuration J belonging to B, by allowing Z to approach J the corre- 
sponding integral can perhaps be taken over J itself, on which v(x, y) is 
continuous. Theorems 1, 5 and 11 below are all proved by this method 
from the formula 


u(x, y) = + f (108 72 ce iat a ") ds, (x, y) in-R, (1) 
v 





after making suitable simplifications. 

THEOREM 1. Let the boundary B of an infinite region R be a finite Jordan 
configuration, and let G(x, y) be Green’s function for R with pole at infinity. 
Then we have 


G(x, y) = SZ logrds + g, OSes I; 


where dc = —dH/2r on B, Jf, do = 1, H(x, y) is conjugate to G(x, y) in R, 
and g is a constant. 
The function 


H(x, y) = J,arg (z — t)do(t), z=x+ v1, 
ts conjugate to G(x, y) in R, locally single valued, and the function 
F(z) = G(x, y) + iH(x, y) = J’, log (2 — t)do(t) + g 


ts analytic in R. 
The critical points of G(x, y) in R are the zeros in R of the analytic function 


F'(z) = | Dias: zin R, 
pez—t 


and thus are the positions of equilibrium in the field of force due to the spread 
o of positive matter on B, where the matter repels with a force proportional to 
the mass and inversely proportional to the distance. 

The last part of Theorem 1 follows by taking the conjugate of F’(z), for 
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the vector 1/(z — #) represents the force at z due to a particle of unit mass 
at f. 

Consideration of the direction of the field of force leads to an analog 
of Lucas’ Theorem: 

THEOREM 2. Under the conditions of Theorem 1, all critical points of 
G(x, y) in R lie in the smallest convex point set Il containing B. No boundary 
point 2 of Tin Rts a critical point of G(x, y) unless B lies on a line through 2. 

The latter part of Theorem 2 is difficult if not impossible to prove by 
previously published methods, as is the latter part of Theorem 3 below. 

2. If in Theorem 1 the boundary B is symmetric in the axis of reals, 
the element of mass do is equal at symmetric points of B. The integral 
expressing F’(z) can be taken over one half of B, using a suitable con- 
vention for any part of B that may lie on the axis, with the integrand 


1 1 
, (2) 
z—t zs-l 








If zis not real and lies exterior to the Jensen circle constructed on the 
segment (f, 2) as diameter, the conjugate of (2) is the force at z due to 
unit particles at ¢ and ?, and whether or not ¢ is real has a non-zero vertical 
component directed away from the axis; we have the first part of 

THEOREM 3. Under the conditions of Theorem 1, let B be symmetric in 
the axis of reals, and let Jensen circles be constructed with diameters the segments 
joining all pairs of symmetric points of B. Then all critical points of G(x, y) 
in R lie in the closed interiors of these Jensen circles. 

A critical point 2 in R not interior to a Jensen circle cannot lie on such a 
circle unless it lies on all Jensen circles for B, that is to say, unless B lies 
wholly on the equilateral hyperbola with vertices 2 and Zp. 

The latter part of Theorem 3 follows from the fact that a point z on the 
Jensen circle for ¢ and ?, the force represented by the conjugate of (2) is 
horizontal. 

At any point 29 = Xo + 7 of R not a critical point, the force in the field 
descrebed in Theorem 1 is orthogonal to the locus G(x, y) = G(xo, yo); for 
if we compute dF/dz, where dz = ¢ dv(\e| = 1) is taken in the direction 
of the normal to this locus in the sense of increasing G, we have dF/dz = 
(0G/dv)/e, OG/Ov > 0, OH /Ov = 0, whence arg [F’(z)] = —arg e, so the 
direction of the force is that of v; the conjugate of F’(z) is the gradient of 
G. This remark obviously applies generally to the field defined by the 
conjugate of an analytic function F’(z), where F = G + 1H. 

By study of the variation in the direction of the force as z traces a level 
locus near B, and also traces a suitably chosen auxiliary Jordan curve in 
R containing K, we study the variation of arg [F’(z)] on these curves, 
and obtain 

THEOREM 4. Under the conditions of Theorem 3 let K be a configuration 
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consisting of a segment a S z SB plus the closed interiors of all Jensen circles 
intersecting that segment. Let a and B lie in R and not be critical points of 
G(x, y). Let K contain k components of B. Then K contains k — 1, k, 
or k + 1 critical points of G in R according as the forces at a and 8 are both 
directed outward, one outward and one inward, or both directed inward. 

3. For forces of the kind described, we prove the 

Lemma. The force at a point P due to a distribution of positive matter 
in a circular region C not containing P is equal to the force at P due to the 
same total mass concentrated at a suitably chosen point of C. 

A circular region is the closed interior or exterior of a circle, or a half- 
plane. 

This Lemma may be proved by inversion in P; the force at P due to a 
particle at Q is in magnitude, direction and sense the vector Q’P, where 
Q’ is the inverse of Q; the proposition is equivalent to the fact that if a 
distribution of positive mass lies in the closed interior of a given circle, 
so does its center of gravity. This Lemma enables us to prove many 
further results, here omitted, from Theorem 1. 

The method of proof of Theorem 1 also yields 

THEOREM 5. Let R be a region of the extended plane whose boundary B 
is finite and consists of disjoint Jordan configurations C and D. Let the 
function u(x, y) be harmonic in R, continuous on R + B, zero on C and unity 
on D. Then in R we have 


u(x, y) = f,logrdo — Jf, logrda + wm, S ds = Sf, de, (3) 


do = —dv/2x on C, do = dv/2x on D, where uy 1s a constant and v(x, y) ts 
conjugate to u(x, y) in R. Wealso have 


v(x, y) = J, arg (2 — t)do — Sf,arg (2 — t)do, 
f(z) =u+t+ w= Slog (z — t)do — J, log (z — t)do + um. 


- The finite critical points of u(x, y) in R are the finite zeros in R of 


FO: ff oy (4) 


Consequently the critical points of u(x, y) in R are precisely the positions 
of equilibrium in the field of force due to a spread o of positive matter on C 
and a spread —o of negative matter on D, where the matter repels with a force 
proportional to the mass and inversely to the distance. 

If R is finite the constant u is zero or unity according as R is separated 
from the point at infinity by C or D; if R is infinite, uw = u(~). 

As an analog of a theorem due to Bécher we prove 

THEOREM 6. Under the conditions of Theorem 5, let C and D lie, re- 
spectively, in disjoint circular regions S and T. Then all critical points of 
u(x, y) in Rliein Sand T. If C and D consist, respectively, of m and n 
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components, then S‘and T contain, respectively, m — 1 and n — 1 critical 
points. 

At a point’ exterior to S and T the force due to the positive mass on 
Cis equal to the force at z due to the same mass concentrated at some point 
Q; of S, by the Lemma, and the force at z due to the negative mass on D is 
equal to the force at z due to the same mass (numerically equal to the total 
positive mass) concentrated at a point Q. of T; the field due to these two 
masses at Q, and Q, has as lines of force the circles through Q; and Qs, so 
the total force at zis not zero. Study of the variation in the direction of 
the force as z traces a circle separating S and T, and traces loci u(x, y) = e 
and 1 — ein R, where ¢ is small but positive, completes the proof. 

4. In Theorems 5 and 6 the constant values zero and unity on Cand D 
may be replaced by any other distinct constant values; such a change 
may modify u(x, y) by an additive constant, which alters um in (3) but 
leaves (4) unchanged, and may modify u(x, y) by a multiplicative constant, 
which multiplies 1%, v(x, y), and hence do by this same constant, and 
leaves (3) and (4) formally unchanged; even if the multiplicative constant 
is negative the conclusion of Theorem 5 persists. This remark will be 
useful in the proof of 

THEOREM 7. Let u(x, y) be harmonic in a region R bounded by a Jordan 
curve J and a Jordan configuration B disjoint from J, continuous on R + 
J + B, zero on J and unity on B. Let the region R, bounded by J and 
containing R be provided with non-euclidean (NE) geometry by mapping 
R, conformally onto the inierior of a circle, and let Il be the smallest NE 
convex set in R, containing B. Then all critical points of u(x, y) in R lie 
in II. 

We make the permissible assumption that J is the unit circle and that 
the origin is a point of R. We extend u(x, y) harmonically across J so 
that u(x, y) is harmonic in the region R’ bounded by B and the reflection 
B’ of Bin J; then u(x, y) is continuous on R’ + B + B’, equal to unity 
on B and minus unity on B’. Both B and B’ are finite. The values of 
do are equal in corresponding points of B and B’, so the integral corre- 
sponding to (4) for R’ can be taken merely over B, with the integrand 


1 1 
s-t z:-—lft 





(5) . 


The conjugate of (5) is the force at z due to unit positive and negative 
masses at ¢ and 1/2, and this force has the direction of the circle through 
z, t and 1/2, in the sense from ¢ toward z. If z is a point of R exterior to 
II, there exists a NE line Z through z disjoint from II, and each vector 
represented by the conjugate of (5) with ¢ on B has a non-zero component 
orthogonal to L directed toward the side of ZL on which II does not lie. 
. Thus z is not a position of equilibrium and Theorem 7 is established. It 
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follows also that no boundary point of II is a critical point unless II reduces 
to a segment of a NE line. 

If B in Theorem 7 is symmetric in a NE line L, there exist precise 
analogs of Theorems 3 and 4, where now the Jensen circles are to be taken 
as NE circles with NE centerson L. Still another related result is 

THEOREM 8. Let u(x, y) be harmonic in a region R bounded by a Jordan 
curve J and two Jordan configurations B, and By, disjoint mutually and from 
J, continuous on R + J + B, + Baz, zero on J, unity on By, minus unity on 
B,. Let the region R, bounded by J and containing R be provided with a 
NE geometry. If there exists a NE line L separating B, and B, in R, then 
the totality of such lines L separates a NE convex set 11, containing B, from 
a NE convex set Il, containing B; all critical points of u(x, y) in R lie in 
Tl, and IIe. 

5. For a region R we denote generically by G(z, 2) Green’s function 
with pole in %, and we study 


= >> AG(z, ax), he >0O.. (6) 
k=1 : 


THEOREM 9. Let R be a Jordan region, and let the points a, a2, ...,; 
am, Of R lie in a smallest NE convex set 11. Then all critical points of u(z) 
defined by (6) also le in I. 

Theorem 9 can be proved by remarking that for large positive M the 
level locus u(z) = M in R consists of m curves, one about each point a;; 
Theorem 7 then applies; as M becomes infinite the smallest NE convex 
set containing the locus approaches II. 

In a similar manner, by considering the level loci u(z) = M and u(z) = 
— M, there may be proved 

THEOREM 10. Let R be a Jordan region, and let aj, ..., my Bi - ++» By 
be points of R which lie in two NE convex sets 1, and Iz separated by every 
NE line separating the a, from the By. Then all critical points of 


u(z) ae MG(z, ax) ~ uiG(z, Br), Aw >O0,ue>0, (7) 


in R lie in Tl, and Mp. 

The situation of Theorem 9 can be studied directly even under a more 
general hypothesis. In Theorem 1 let 2% be a finite point of R; prises 
u(x, y) in (3) as G(z, 2), with C as B and D as the locus G(z, %) = 
When WM varies, v(x, y) does not ay. and we have {, do = 1, so as _ 
becomes infinite we obtain 


G(z, 2) = J, log r do — log|z — | + wo, S do = 1, 
F(z) = G(z, 2%) + iH(z, 2) = J log (z — t)do — log (z — 2) + um, 
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THEOREM 11. Let R be a region bounded by a finite Jordan configuration 
B, and let ay, ..., am be finite points of R. Then the critical points of u(z) 
defined by (6) are the positions of equilibrium in a field of force due to a spread 
o of positive matter of total mass =), on B, and to negative particles of re- 
spective masses— dx at the az. 

One of many consequences of Theorem 11 is 

THEOREM 12. Under the conditions of Theorem 11 let B lie in a circular 
region S and the a, in a circular region T disjoint from S. Then all critical 
points of u(z) in R lie in S and T. 

6. As is well known, equation (1) expresses u(x, y) as the sum of a 
simple potential and a double layer potential; in any case we may write 
from (1) at least formally? 


—1 1 
u(x, ¥) = 2S) log r dv + ee (2 — #)] 
Qe Js Qn J » 
5 1 1 
= Ff te + do + 5 [wang — Ola — 5 f are (@ — du, 
28 SF Qa Qa 


wt io = > ff tog (@ — fav — 5 [log (@ — Da + 
2r J p 2r 


a 

if log (2 — é)du, 

2a J 2 

do $|. & 4 du 
f{4-dssh+ess 8) 
p2—t whe—tls ws,2—!t 

The gradient of u(x, y) is the conjugate of f’(z), for which there is an 
obvious interpretation as a field of force due to spreads o and iu of matter 
over B and particles at the end-points of B. Theorems 1, 5 and 11 reduce 
this general configuration to simpler manageable situations for suitable 
u(x, y), and those theorems by no means exhaust the possibilities. For 


instance, if u(x, y) has the constant value unity on an are A: (a, 8) of B, 
the corresponding last two terms of (8) reduce to 


oa. ene See 

—_ | ee ; (9) 

2r lz —B 2-a 
the conjugate of (9) is a field of force due to skew particles at B and a, of 
masses +1 and —1, respectively; these particles exert forces at z equal 
to the inverse distance, in the directions arg ( — 8) + m/2 and arg (z — 


a) — 4/2. The lines of force in this field are the circles of the coaxial 
family determined by a and @ as null circles. Such fields of force may be 





F(z) 











f'(2) 
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combined with similar ones, or with such fields as occur in Theorems 1, 
5 and 11. 

Moreover, if u(x, y) is harmonic interior to a circle C and is extended 
harmonically across an arc of C on which u(x, y) vanishes, 0u/Ov has equal 
and opposite values on the two sides of the complementary arc of C, so 
over the latter arc the first part of the integral in (1) or (8) vanishes. 
Several results are consequences of these remarks: 

THEOREM 13. Let the points ay, Bi, a2, Be, ..., On = a0, Bn = Bo lie on 
the unit circle C: 2 = e” in the counterclockwise order indicated, and let the 
function u(z) harmonic and bounded interior to C take a positive constant 
value y;. on each arc a, << 0< B, and the value zero on each arc By < O< ati. 
Let NE lines Ao, Ai, ..., An-1 be so chosen that A, is an.arc of a circle be- 
longing both to the coaxial family determined by a, and B, and to that determined 
by a+, and By+1. Then the NE convex set interior to C bounded by the A; 
contains all critical points of u(z). 

These arcs A; are special cases of W-curves, which are of importance in 


the entire study of critical points, namely, the locus of points at which the 


forces due to two groups of particles have the same direction but opposite 
senses; the groups of particles (here the » groups a;, 8x) have relations of 
equality of mass between the elements of a group, but between groups all 
possible relative masses of constant sign are considered. As in Theorem 
13, the W-curves form the boundaries of regions necessarily containing the 
critical points. A double layer distribution on an arbitrary curve C is 
the limit of a double layer distribution piecewise constant on C as in 
Theorem 13; the lines of force for the original distribution are the circles 
normal to C, and the W-curves are the circles normal to C in two points 
plus the locus of points 2 from which two circles tangent at 2) and normal 
to C can be drawn. 

THEOREM 14. Let R be a region bounded by the unit circle C and by a 
Jordan configuration B interior to C. Let the function u(x, y) be harmonic 
and bounded in R, continuous on R + B + C except in the end-points of an 
arc A of C, unity on B and the interior points of A, zero on the interior points 
of C—A. Let Il bea NE convex set containing B, and let the set 1, consist 
of II plus all points of arcs joining A and Il, of circles of the coaxial family 
determined by the end-points of A as null circles. Then all critical points of 
u(x, y) in R lie in Th. 

The conclusion of Theorem 14 applies also to an arbitrary linear combina- 
tion with positive coefficients of the two harmonic measures w (z, A, | z | < 1) 
and w(z, B, R). 

Theorem 14 is a sharpened analog of Theorem 7; a similar sharpened 
analog of Theorem 8 also exists, and as in Theorem 13 we may replace a 
single arc of C by a sum of arcs of C; both positive and negative values of 
u(x, y) on B and C may be admitted. Indeed, all the results formulated 
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in the present note are intended only as illustrations of an extensive theory 
which the writer plans to develop elsewhere. 


1 For detailed references to the literature, the reader may consult a forthcoming 
book: Marden, M., Geometry of the Zeros of a Polynomial, to be published as a volume 
of Mathematical Surveys by the American Mathematical Society. 


2 Equation (8) is Cauchy’s integral formula for f’(z), involving an unusual term to 
include discontinuous u(x, y). 


SOME PROPERTIES OF ROTATIONAL FLOW OF A PERFECT 
GAS 


By P. NEMENYI AND R. PRIM 
NAVAL ORDNANCE LABORATORY, WHITE OAK, MARYLAND 


Communicated by S. Lefschetz, January 14, 1948 


1. Introduction—Complete analytic solutions of definite boundary 
value problems of compressible rotational flow are extremely difficult to 
obtain. Therefore, in the present stage of study of rotational gas flow the 
inverse, or semi-inverse, approach promises to prove fruitful. We propose 
to deduce and to study flow patterns satisfying the differential equations 
of gas flow and having throughout the field certain prescribed geometric, 
kinematic or physical properties; but will not make them obey any pre- 
scribed boundary conditions. 

Usually entire systems of flow patterns, rather than individual solutions, 
result from this type of approach. If an imposed condition is proved not 
to be satisfied by any possible flow pattern, then this negative result is of 
significance, as it establishes a ot property of all flows satisfying the 
equations of the problem. 

In various fields of fluid mechanics the semi-inverse approach has been 
used with considerable success. It is sufficient to recall the contributions 
of Beltrami, Massotti, Jeffreys, Hamel, Oseen, Kampé de Fériet,' Bateman 
and Tollmien.? None of these investigations deals however with the rota- 
tional flow of a perfect gas. : 

The reason for this is probably that such an investigation—aunless re- 
stricted to cases of uniform stagnation enthalpy—requires complicated 
and difficult eliminations if based on the familiar equations of the problem. 

In the recent investigations of Munk and Prim* ‘ this elimination is 
accomplished with complete generality as far as steady flow of a perfect 
gas is concerned. This makes the application of the semi-inverse ap- 
proach to rotational gas flow practicable. Therefore the results of Munk 
and Prim are used throughout the present paper. The present paper is 
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a general preliminary report on our investigations, parts of which will be 
published in detail elsewhere. 

2. Some Properties of the Plane Flow of a Perfect’ Gas.—A natural 
coérdinate system for the analysis of problems of plane fluid motion employs 
as codrdinate curves the streamlines (¢ = const) and their orthogonal 
trajectories (7 = const). In this codrdinate system the vorticity of the 
reduced velocity field is expressed by 


‘ 


|curl @| hae w (1) 
curl @| = —— —g 
Bigs OE 
(g: dé and godn are components of the vector element of arc length). 
The vorticity equation and the continuity equation‘ take on the following 





form: 
re) w? On (gow) ieee " 
area dé oe. @) 
= [g.(1 — w?)’w] = 0. (3) 
On 


These equations enable us to settle certain general questions concerning 
plane gas flows. 

We ask first under what conditions the velocity is constant along each 
streamline. Since the ultimate velocity is in every case constant along 
each streamline, our question is equivalent to asking under what conditions 
the reduced velocity is constant along each streamline. 

Introduction of the requirement w = w() inte equations (2) and (3) 
yields the geometric conditions: 


07 In go 

is 0 or ge = a(€)b(n) 

o 1 

ee we. ae = ht) 
On 


An analysis of this geometric limitation leads to the following theorem: 

THEOREM 1. In plane steady flow of a perfect gas, in absence of an external 
field of mass force, the velocity magnitude can have a constant value along each 
streamline only if the streamlines are concentric circles or parallel straight lines. 

Obviously, in these cases in which w = const (and hence also v = const) 
along each streamline, curl ®@ and curl 7 also have a constant magnitude 
along the same lines. The question arises whether the converse theorem 
is also true, that is, whether the only two dimensional flow patterns possible 
in a perfect gas such that the streamlines are-lines of equal vorticity are 
those in which they are also lines of equal velocity magnitude. 
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To settle this question we consider first flow with constant reduced 
vorticity magnitude along each streamline. Crocco’s pressure theorem as 
reformulated® for a general]- distribution of stagnation enthalpy requires 
that the pressure remain constant along each streamline if the reduced 
vorticity does so; that is, p = p(é). On the other hand, according to a 
familiar relation ; 


p 
po(E) 





me (1 iene w?)b+1 


(po denotes the stagnation pressure, which remains constant along each 
streamline between shock fronts). Hence ‘curl D| = f(&) implies that 
w must be also a function of & alone. 

A very similar argument shows that if the lines of equal velocity magni- 
tude (isovels) and the lines of equal vorticity (isocurls) of the reduced 
velocity field coincide, then they coincide also with the streamlines. 

Starting from the vector identity 


curl = a(é) curl ® + grad a(t) X @ 


we have transformed the above statements relating to the reduced velocity 
field into analogous statements concerning the actual velocity field and 
proved the following theorems: 

THEOREM 2. In plane steady flow of a perfect gas, in absence of an external 
field of mass forces, constant vorticity along each streamline implies constant 
velocity along the same lines. 

THEOREM 3. In plane steady flow of a perfect gas, in absence of an external 
field of mass forces, the coincidence of the isovels and the isocurls implies that 
they coincide also with the streamlines. 

Our three theorems prove that all flow fields having coincident isovels 
and streamlines, or coincident isocurls and streamlines, or coincident 
isovels and isocurls are restricted to flows the streamlines of which are 
concentric circles or parallel straight lines. 

In a paper by R. Prim® the question is investigated whether the stream- 
line pattern can form, with its orthogonal trajectories, an isothermal net, 
that is a net in which g; = go. The following theorem is established. 

THEOREM 4. In plane steady flow of a perfect gas in absence of an external 
field of mass forces, the streamlines and their orthogonal trajectories can form 
an tsothermal net only if the streamlines are concentric circles, parallel straight 
lines or radial straight lines: the last case being possible only if the reduced 
velocity field is irrotational. 

It will be instructive to compare the geometric restrictions imposed by 
the above four conditions upon the pattern of gas flow with the corre- 
sponding restrictions imposed upon incompressible non-viscous and in- 
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compressible viscous flow. 
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3. On Flow Fields for Which, in a Cylindrical Isothermal Codérdinate 


System All Three Velocity Components Depend Only on One Codrdinate.— 
For irrotational, plane gas flow this question was investigated by- Toll- 
mien.2 We have extended the study to rotational flow, at the same time 
allowing an additional velocity component parallel to the elements of the 
cylindric eoérdinate surfaces. In one respect, however, our investigation 
is less general than that of Tollmien: we have investigated only systems 
of such cylindrical coérdinate surfaces which in the xy plane form an 
isothermal net, while Tollmien admitted orthogonal nets of any kind. 

The main result of the investigation can be summarized in the theorem: 

THEOREM 5. The only isothermal cylindrical system of coordinates for 


which the reduced velocity components of any steady non-parallel flow of a 
perfect gas, in absence of an external field of mass forces, are functions of one 
of the isothermal coordinates only are farmed by logarithmic spirals and their 
Thus the coérdinate systems satisfying our kinematical 


limiting cases. 
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conditions are exactly the same as those found by Tollmien for plane 
potential flows. With reference to this spiral-cylindric coérdinate system 
the differential equations can be reduced to a pair of ordinary non-linear 
differential equations of the first order. 

4. Massotti Flow.—Massotti sought incompressible, non-viscous flow 
patterns for which v, = 0, and for which at the same time (curl 7), = 0. 

There is a wide variety of compressible flow patterns the reduced velocity 
vector of which satisfies the corresponding relations w, = 0 and (curl ®), = 
0. We satisfy the latter equation by setting 


_ 2, op 
Ox’ ~ Oy 


“-(5) Ge 


and the continuity equation takes on the form: 


, Od Qw? dp dw? 
Arw[1 — w*] = B| — — — |, 
ol ed : ks Qe Dy =| 


Wz 


Hence 





® X curl @ 

1 — w? 
as ® X curl @ is obviously perpendicular to the x, y plane and equal to 
*/2(Ow*/dz). 

Hence the only condition imposed upon ¢ is the continuity equation 
which is identical with that of irrotational plane compressible flow. We 
can immediately see that a wide variety of Massotti flows can be built from 
families of plane irrotational flows in the following way. Assume (x, y, 
C1, C2, ... €n) is the potential of the reduced field of an n-parameter family 
of irrotational flows. Replace each of these parameters by an arbitrary 
function of z; the ¢ thus obtained can be considered as the potential from 
which the x and y components of families of reduced velocity fields of the 
Massotti type can be derived. 

For example, we have based upon the simple irrotational Prandtl-Meyer 
flow a pseudo-plane rotational flow past a developable surface. 

5. Beltrami Flows and Other Investigations —Beltrami’ was the first 
to investigate velocity fields which satisfy the condition » X curl v = 0. 
In recent memoranda® * ' the present writers investigate flows of perfect 
gases which satisfy this relation (Beltrami Flow Proper). They found it 
fruitful to study also the broader family of flows for which ® X curl ® = 0 
(Generalized Beltrami Flow). 

Further semi-inverse investigations under way include the study of gas 


Our vorticity equation, curl ] = 0, is automatically satisfied 
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flow patterns for which all streamlines are straight lines or for which the 
hodograph space degenerates into a surface or a curve. In addition, the 
basis of these semi-inverse methods, the equations of Munk and Prim, are 
being generalized to include certain classes of non-steady flow and of com- 
pressible substances with a more general state equation. 


1 See, for example, Berker, A. Ratib, “Sur quelques cas d’integration des équation du 
mouvement d’un fluide visqueux incompressible,’ Paris et Lille, 1936. 

2See Tollmien, W., “Zum Ubergang von Unterschall—in Uberschallstrémungen,” 
Z. f. Angewandte Mathematik und Mechanik, 17, 117-136 (1937). 

3 Munk, M. M., and Prim, R., “‘On the Multiplicity of Steady Gas Flows Having 
the Same Streamline Pattern,”’ Proc. Nat. Acad. Sci., 33, 187-141 (1947). 

4 Munk, M. M., and Prim, R., “‘On the Canonical Form of the Equations of Steady 
Motion of a Perfect Gas,’”’ Naval Ordnance Laboratory Memorandum No. 9169 (June, 
1947). 

5 Prim, R., ‘‘Extension of Crocco’s Theorems to Flows Having Non-Uniform Stagna- 
tion Enthalpy,’’ Naval Ordnance Laboratory Memorandum No. 9286 (August, 1947); 
Phys. Rev., Jan. 15, 1948. 

6 Prim, R., ‘‘On the Existence of Steady Gas Flow in Plane Isothermal Streamline 
Patterns,’’ Naval Ordnance Laboratory Memorandum No. 9267 (August, 1947); Bull. 
Amer. Math. Soc. (forthegming). 

7 Beltrami, E., ‘“‘Considerazioni idrodinamiche,’’ Rendiconti Istituto Lombardo, 22, 
121-130 (1889). 

8’ Neményi, P., and Prim, R., ‘“Some Patterns of Vorticose Flow of a Perfect Gas,” 
Naval Ordnance Laboratory Memorandum No. 9219 (June, 1947). 

® Neményi, P., and Prim, R., ‘‘Some General Properties of Beltrami Flows of a 
Perfect Gas,’’ Naval Ordnance Laboratory Memorandum No. 9397 (October, 1947). 

10 Neményi, P., and Prim, R., “‘Some Interesting Cases of Steady Helicoidal Beltrami 
Flow of a Perfect Gas,’’ Naval Ordnance Laboratory Memorandum No. 9426 (No- 
vember, 1947). 


REPRODUCTIVE DIAPAUSE IN DROSOPHILA ROBUSTA 


By Hampton L. CARSON AND HARRISON D. STALKER 
DEPARTMENT OF ZOOLOGY, WASHINGTON UNIVERSITY, St. Louts, Mo. 
Communicated by Th. Dobzhansky, January 19, 1948 


Introduction.—Studies of microevolutionary processes in natural popu- 
lations of Drosophila have been hampered by a lack of ecological data. For 
most species of the genus, the basic facts concerning population size, mi- 
gration, longevity, specific breeding sites, length of breeding season and 
nature of overwintering populations are unknown. Workers in this field 
have not ignored these subjects * *, but it is clear that preoccupation with 
other basic problems has generally held up experiments specifically planned 
to answer these questions. It is, nevertheless, true that a full understand- 
ing of natural selection, genetic drift and other processes, as they are oper- 
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ative in natural populations of Drosophila, will not be attained until the 
ecology of the various species concerned is better known. The data pre- 
sented in this paper are based largely on observations which were incidental 
to a study of the gene arrangements and morphology of Drosophila robusta, 
throughout the warm seasons of the year in a single locality near St. Louis, 
Missouri. 

Although of a preliminary nature, these data show that a striking repro- 
ductive arrest, or diapause, sets in during the late summer and early au- 
tumn in the natural population, despite the fact that temperatures ap- 
parently favorable for the breeding of the flies in nature continue ordinarily 
until the middle of October. The importance of this fact for the interpre- 
tation of population dynamics in this species is obviously great, and further, 
more detailed investigations of this and similar subjects are planned. 

Observations.—The specimens of Drosophila robusta on which the following 
observations were made were collected from paper cups baited with ripe 
banana mash and set in a particular deciduous woods, supporting consider- 
able undergrowth, near Olivette, St. Louis County, Missouri. In 1946, 
flies were collected over nine periods of variable length between March 5 
and December 10 (see first column of table 1). As the offspring of these 


TABLE 1 
FERTILITY OF MATURE Wip-Caucut FEema.es oF Drosophila robusta. 1946 COLLECTION 


FEMALES PRODUCING OFFSPRING 


No. OF MATURE* WITHIN SIX DAYS 

COLLECTION PERIOD FEMALES COLLECTED No. PER CENT 
I. Mar. 5-Apr.6 ° 38 27 71.1 
II. Apr. 18—May 15 99 69 69.7 
III. May 18-Jun. 1 153 120 78.4 
IV. Jun. 17-Jun. 19 ° 203 171 84.2 
V. Jul. 3-Jul. 9 108 84 77.8 
VI. Jul. 23—Aug. 6 75 64 85.3 
VII. Aug. 15-Sept. 3 42 29 69.0 
VIII. Sept. 13-Sept. 24 52 2 3.9 
IX. Oct. 14-Dec. 10 45 0 0.0 


* Wild-caught females which had obviously only recently emerged from the pupal 
case are not included in this table. 


wild flies were to be used for the determination of the frequencies of the 
various gene arrangements and an F2 generation from the majority of the 
wild flies obtained, the following procedure was adopted. Each wild- 
caught female was placed individually into a separate culture vial soon after 
capture; after three days it was transferred to a fresh culture vial and at 
the end of six days both vials were examined for the presence of larvae. If 
larvae were present, the female was recorded as having been inseminated 
in nature and the larvae were raised, some being used for the cytological 
determination of the salivary gland chromosome sequences and the re- 
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mainder saved as a mass F; from the original wild female. If no larvae 
were observed in six days, the female was mated to a male captured on or 
about the same date and the F; handled in a similar manner. From the 
middle of June on, a considerable number of females were captured which, 
because of their light coloration, had obviously only recently emerged from 
the pupal case. These were separated from the ‘“‘mature’’ ones, aged for 
ten days en masse and then mated individually to males captured on the 
same day. 

Table 1 shows the total number of females, judged as relatively mature, 
which were caught during the nine collection periods and the per cent of 
these which were found to have been inseminated in nature, as indicated by 
the production of offspring in the laboratory within six days. It will be 
seen that this percentage remains high from March 5 through August 6, 
drops in the next period (August 15 through September 3) and falls abruptly 
to a very low figure for females collected after September 13. It is prob- 
able, from the 1947 data given below, that most of the flies failing to pro- 
duce offspring within six days were virgin, although from the 1946 data 
alone, the possibility of some sort of slowing-down of egg production by in- 
seminated females cannot be excluded. With a few exceptions, however, 
all the females collected on or after September 13 produced offspring after 
being mated to males captured at the same time, although in many cases the 
production of fertile eggs by these females was delayed for a week or more 
after being placed with a male. 


TABLE 2 


ELAPSED TIME BETWEEN EMERGENCE AND THE FIRST APPEARANCE OF F2 LARVAE FOR 
F, Fires From Witp-Caucnt Drosophila robusta. 1946 COLLECTION 


NO. OF DAYS ELAPSING FROM EMERGENCE OF 


COLLECTION PERIOD OF No. OF Fis F\ FLIES UNTIL F2 LARVAE WERE PRODUCED 
P\ WILD FLIES OBTAINED MEAN MODE RANGE 
IV. Jun. 17—Jun. 19 146 14.5 13 7-31 
V. Jul. 3-Jul. 9 150 17.2 14 8-50 
VI. Jul. 23—Aug. 6 106 19.5 15 7-51 
VII. Aug. 15-Sept. 3 94 25.4 21 10-56 
VIII. Sept. 13—Sept. 24 109 33.5 27 11-65 
IX. Oct. 14—Dec. 10 31 33.0 34 13-55 


F, generations were obtained in the following manner from most of the 
above wild-caught females or pairs of wild flies. F; flies of both sexes were 
maintained in mass culture in vials and changed every three days until F, 
larvae were observed. For each of these Fis, the date of the latter event, as 
well as that of the emergence of F; flies was recorded as a matter of routine. 
These data are thus useful for an analysis of the speed with which the F; 
flies reached a condition of active reproduction as indicated by the pro- 
duction of offspring. Table 2 shows that the F,s from wild flies caught 
from August 15 on attained this condition much more slowly than those 
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collected in June, July and the first week in August, Correspondence be- 
tween these data and those for the fertility of the wild-caught parental flies 
(table 1) is striking, and it is clear that the reproductive delay manifested 
by the parental flies caught on or after August 15 is reflected in their off- 
spring, despite the fact that the latter were kept under uniform-laboratory 
conditions of temperature (25.5°C.). 

The above data appear to be significant in view of the fact that in old 
laboratory stocks of Drosophila robusta, sexual maturity (as indicated by 
willingness to mate) is attained about four days after emergence for females 
and about eight days for males. In this species, mating appears to occur 
several days before the females produce offspring, although further study 
of this characteristic in laboratory stocks is desirable. 


TABLE 3 


INSEMINATION OF WILD FEMALE Drosophila robusta CAPTURED DURING SEPTEMBER 
AND OCTOBER, 1947 


NO. OF FEMALES 


TOTAL NO. OF INSEMINATED IN 
DATE OF CAPTURE ° FEMALES NATURE 
1. Insemination determined by observation of hatching eggs: 
Sept. 28—Oct. 2 14 0 
Oct. 3-5 30 8 
Oct. 7-9 27 1 
Oct. 10-12 49 1 
Oct. 138-15 28 0 
2. Insemination determined by dissection: 
Oct. 14-15 18 1 
Oct. 16-18 25 1 
Oct. 20-22 21 1 
Oct. 23-24 : 17 0 
Oct. 25-27 26 0 
Total ‘ 256 13 
Per cent inseminated in nature 5.1 


In 1947, collections were made during two periods from the same locality 
in the same woods. These periods extended from May 19 through June 7 
and from September 28 through October 27, respectively. In order to re- 
duce the labor involved and to expedite analysis of the chromosomes, each 
wild-caught female in the spring collection was placed with a male immedi- 
ately upon capture. For this reason, the spring data are of no use in the 
analysis of fertility in nature. All but very few pairs produced offspring in 
six days; no Fy, generations were raised in 1947. 

The fall (1947) collection was analyzed as follows. Females were placed 
individually into culture vials shortly after capture. At three-day inter- 
vals, each was changed to a fresh vial and the surface of the food closely 
examined under the binocular microscope for the presence of eggs or larvae. 
The date of the appearance of eggs was recorded, and if no larvae hatched 
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from these in the six days following, the female was considered to be virgin 
and was mated to a wild-caught male. In addition to these, 107 females 
caught between October 14 and October 27 (when the collecting season was 
ended by the abrupt onset of cold weather) were dissected immediately on 
being broufght to the laboratory. For each of these, notes were made on the 
size of the ovary, the character and amount of body fat and whether the fe- 
male had been inseminated in nature, as indicated by the presence of sperm 
in the ventral receptacles or spermathecae. 

Table 3 shows that approximately 95% of the females caught during this 
period were virgin, as determined either by watching the fate of eggs laid 
in the laboratory (upper part of table) or by dissection (lower part of table). 
Of the 107 females which were dissected, furthermore, none had functional 
ovaries as indicated by the presence of large filamented ovarian eggs, 
although in two of them the ovaries were medium in size and contained 
filamented eggs approximately two-thirds mature size. In almost all the 
females, the ovaries were very small and rounded (approximately 380 u in 
length), similar to those of newly emerged flies. _A high percentage of the 
flies, furthermore, had heavy sheets of body fat, sometimes almost com- 
pletely filling the abdomen. Such a condition is quite unlike that in a 
newly emerged fly in which the residual pupal fat, in the form of spherical 
“nodules” is easily recognized. In 13 of the females, nodular fat was 
found; this was taken as indicative of relative youth, although, in some of 
these, fat was also present in small sheets. 

Of particular interest were the three dissected flies which had been in- 
seminated in nature. Like most of the virgins, these flies had abundant 
body fat present in sheets and very small ovaries. The latter, however, 
were peculiarly elongated and flattened apically, while proximally they pre- 
sented a slightly swollen appearance. This swollen area contained small 
ovarian ova of a distinct yellowish-brown color, quite unlike the colorless or 
whitish ova usually encountered. Ovaries of this unusual type were also 
found in one of the virgin flies. ; 

Of the 148 females collected between September 28 and October 15, 10 
had been inseminated in nature and 138 were virgin. The elapsed time be- 
tween capture and the laying of the first infertile eggs by the virgins was 
recorded. This was extremely variable, ranging from 2—25 days, the mean 
number of days elapsing being 10.3. The 10 inseminated females, further- 
more, required a similar mean number of days (9.1) before fertile eggs were 
laid. One of these females showed an elapsed time of 25 days from the time 
of capture until the first eggs were laid. 

Conclusions—The 1946 data show clearly that, beginning as early as 
August, wild females of Drosophila robusta captured in the fall reproduce 
much more slowly than those collected in spring and early summer... A 
similar delay is manifested in the offspring of these flies. The 1947 data 
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support this conclusion and suggest that the fall population of females, at 
least, consists principally of virgin flies with undeveloped ovaries and large 
quantities of body fat. At least some of the inseminated females have very 
small, peculiar ovaries which suggest regression from a former more highly 
developed and possibly functional condition. This relationship between 
the size of the ovary and the amount of body fat closely resembles the con- 
dition in the reproductive diapause which has been shown to occur in other 
insects in the fall (e.g., ‘dissociation gono-trophique”’ in Anopheles maculi- 
pennis*). 

A somewhat similar diapause has been described for Drosophila nitens.5 
Laboratory strains of this species, despite the fact that they were kept at 
25°C., suspended reproduction from November to March for two consecu- 
tive years. Bertani® has shown that this reproductive diapause can be 
broken by subjecting the adult flies to temperatures from 2-5°C. for 10 
days. In Drosophila robusta, on the other hand, no interruption or delay of 
breeding in old laboratory stocks has been observed and a mean of ten days 
of laboratory conditions appears to be sufficient to break the diapause in 
wild-caught females. 

Biologically, the inference to be drawn is that Drosophila robusta over- 
winters as an adult and that this hibernation period is preceded by a physio- 
logical change-over from egg production to the deposition of body fat. 
The factors, whether genetic or environmental or both, which may be re- 
sponsible for the initiation of the diapause are unknown and will be the ob- 
ject of further investigation. From the point of view of population dy- 
namics, however, the centtfal fact which emerges from these observations is 
that the fall population of this species consists largely of non-breeding 
individuals. 


1 Dobzhansky, Th., and Epling, C., Carnegie Inst. Wash. Publ. No. 554, pp. 1-46 
(1944). 

2 Dobzhansky, Th., and Wright, S., Genetics, 28, 304-340 (1943). 

3 Spencer, W. P., Ohio Jour. Sci., 41, 190-200 (1941). 

4 Swellengrebel, N. H., Ann. Inst. Pasteur, 43, 1370-1389 (1929). 

5 Buzzati-Traverso, A., Istituto Lombardo di Science e Lettere, Rendiconti, 77, 37-49 
(1944). 

6 Bertani, G., Nature, 159, 309 (19-47). 
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